AD-A144  445  JOINT  SERVICES  ELECTRONICS  PROGRAPKU)  STANFORD  UNlV  Cfl 
STANFORD  ELECTRONICS  LABS  S  F  LUNDSTROM  JUN  84 
DAAG29-81-K-0057 


UNCLASSIFIED 


F/G  975 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS -1963-A 


1C  FILE  COPY  AD-A144  445 


STANFORD  ELECTRONICS  LABORATORIES 

I  . .  ~ 

OERARTIKNT  OF  ELECTRICAL  ENGINEERING 

STANFORD  UNIVERSITY  •  STANFORD.  CA  94305 


VNNUAL  PROGRESS  REPORT  NO. 


April  1983  through  31  March  1984 


This  work  was  supported  wholly  by  the 
Joint  Services  Electronics  Program 
(U.S.  Army,  U.S.  Navy,  and  U.S.  Air  Force) 
Contract  DAAG-29-81-K-0057 


Reproduction  in  whole  or  in  part 
is  permitted  for  any  purpose  of 
the  United  States  Government. 


This  document  has  been  approved 
for  public  release  and  sale; 
its  distribution  is  unlimited. 


84  08  07  106 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Data  Entargd! 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  NUMBER 


4.  TITLE  (and  Subtitlal 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


RECIPIENT  S  CATALOG  NUMBER 


S.  TYPE  OF  REPORT  &  PERIOO  COVERED 


Joint  Services  Electronics  Program  Annual  Report 


7.  AUTHORU) 


Annual  Report 
April  1983-March  1984 


8.  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NUMBERU) 


Dr.  Stephen  F.  Lundstrom,  Program  Director. 


9.  PERFORMING  ORGANIZATION  NAME  ANO  ADDRESS 

Stanford  University 

Stanford  Electronics  Laboratories 

Stanford,  California  94305 


11.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

U.S.  Army  Research  Office 
P.  0.  Box. 12211 

Research  Triangle  Park,  NC  27709 


14.  MONITORING  agency  NAME  &  ADDRESS  III  dIH.  from  Controlling  Offieal 

Same  as  11  plus  for  some  duties: 

Office  of  Naval  Research,  Stanford  Branch-Durand  165 
Stanford  University,  Stanford,  CA  94305 


16.  DISTRIBUTION  STATEMENT  (of  this  report! 


DAAG29-81-K-0057 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  &  WORK  UNIT  NUMBERS 


12.  REPORT  DATE  13.  NO.  OF  PAGES 

June  1984  84 


IS.  SECURITY  CLASS.  (o«  thi«  report) 


Unclassified 


ISa.  DECLASSIFICATION /DOWNGRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  abftract  entered  in  Block  20.  if  different  from  report! 


is.  supplementary  notes  views  opinions,  and/or  findings  contained  in  this  report 

are  those  of  the  author(s)  and  should  not  be  construed  as  an  official  Department  of 
the  Army  position,  policy,  or  decision,  unless  so  designated  by  other  documentation, 


19.  KEY  WORDS  (Continue  on  reverse  tide  if  necetfery  end  identify  by  block  number) 


20.  ABSTRACT  (Continue  on  reverse  tide  if  necetfery  end  identify  by  block  number) 

This  is  the  third  annual  report  of  the  research  under  JSEP  contract 
No.  DAAG-81-K-0057 . 


DD,r,1473 

EDITION  OF  I  NOV  BE  IS  OBSOLETE 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Whan  Oata  Entaradl 


Stanford  University 


Stanford  Electronics  Laboratories 


Table  of  Contents 


1.  VERY  LARGE  SCALE  INTEGRATION  6 

1.1  Properties  of  Materials:  Application  of  Channeling  Radiation  to  a  Study  of  the  Properties  6 

of  Materials 

1.1.1  Scientific  Objectives  6 

1.1.2  Progress  7 

1 .1 .2.1  Channeling  Radiation  from  Si  with  an  Oxygen  Platelet  impurity  7 

1 . 1 .2.2  Calculated  Potentials  and  Eigenvalues  in  GaAs,  AIGaAs,  and  Al As  11 

1 . 1 .2.3  Low  Temperature  Channeling  Radiation  1 4 

1.1 .2.4  Electron  Channeling  Radiation  from  LiH  and  LiO  14 

1. 1.2.5  12.6  MeV  Electron  Channeling  Radiation  from  Diamond,  Si  and  LiF  17 

1.1.3  References  21 

1.1.4  Recent  JSEP  Publications  21 

1.2  PHYSICS  AND  TECHNOLOGY  OF  SUBMICRON  DEVICES 


1.2.1  Electrical  Performance  and  Physics  of  isolation  Region  Structures  for  VLSI 

1.2. 1.1  Scientific  Objectives 

1.2. 1.2  Progress 

1 .2.1.3  Summary  of  Results 

1.2. 1.4  References 

1 .2. 1 .5  Recent  JSEP  Publications 

1.2.2  Alternative  Gate  Dielectrics  for  Submicron  MOS  VLSI 

1. 2.2.1  Scientific  Objective 

1. 2.2.2  Progress 

1. 2.2.3  Summary 

1. 2.2.4  References 

1. 2.2.5  Recent  JSEP  Publications 

1.2.3  Electronic  Properties  of  Polycrystalline  Silicon 

1. 2.3.1  Scientific  Objective 

1. 2.3.2  Progress 

1. 2.3.3  Summary 

1. 2.3.4  References 

1. 2.3.5  Recent  JSEP  Publications 

1.2.4  Ohmic  Semiconductor  Contacts  with  Low  Minority-Carrier  Recombination 
Velocity 

1. 2.4.1  Scientific  Objectives 

1. 2.4.2  Progress 

1 .2.4.3  Future  Plans 

1. 2.4.4  Recent  JSEP  Publications 

1.2.5  Interaction  of  Transition  Metals  with  Silicon  Surfaces 

1. 2.5.1  Scientific  Objective 

1. 2.5.2  Progress 

1. 2.5.3  References 

1. 2.5.4  Recent  JSEP  Publications 

1.2.6  Interaction  of  Impurities  with  Silicon  Surfaces  and  interfaces 

1. 2.6.1  Scientific  Objective 

1. 2.6.2  Progress 

1 .2.6.3  References  67 


JSEP  Annual  Report 


Page  i 


June  1984 


888888888288  8822gg88S388fc&£>fc&8!228 


Stanford  University 


Stanford  Electronics  Laboratories 


1.2.7  Interconnections  for  High  Density,  High  Performance  VLSI 

1. 2.7.1  Scientific  Objectives 

1. 2.7.2  Progress 

1.2.7 .3  References 

1.2.8  Study  of  Compound  Semiconductor  Device  and  Materials  Physics 

1. 2.8.1  Scientific  Objective 

1. 2.8.2  Progress 

2.  INFORMATION  SYSTEMS 

2.1  Real  Time  Statistical  Signal  Processing 

2. 1 . 1  Scientific  Objectives 

2.1.2  Progress 

2.1. 2.1  Determining  the  Number  of  Sources 

2.1 .2.2  Coherent  Signals  and  Interference 

2.1 .2.3  Other  Applications  of  Eigenstructure  Methods 

2.1 .3  Summary  of  Results 

2.1 .4  References 

2.1 .5  Recent  JSEP  Publications 

2.2  Signal  Processing  and  Compression 

2.2.1  Scientific  Objectives 

2.2.2  Progress 

2.2.2. 1  Vector  Quantization 

2.2.2.2  Advances  in  Cryptography 

2.2.2.3  Advances  in  Data  Compression  and  Information  Inequalities 

2.2.2.4  Data  Compression  for  Computer  DeteStructures 

2.2.3  References 

2.2.4  Recent  JSEP  Publications 


JSEP  Annual  Report 


Pageii 


June  1984 


288222j2iSJ33i5!3i2222oi3  «  2  2  2  8  8  3  9 


Stanford  University 


Stanford  Electronics  Laboratories 


List  of  Figures 


Figure  1-1: 
Figure  1-2: 

Figure  1-3: 

Figure  1-4: 

Figure  1-5: 

Figure  1-6: 
Figure  1*7: 
Figure  1-8: 
Figure  1-9: 
Figure  1-10: 

Figure  1-11: 

Figure  1-12: 

Figure  1-13: 

Figure  1-14: 

Figure  1-15: 

Figure  1-16: 

Figure  1-17: 

Figure  1-18: 
Figure  1-19: 


High  resolution  TEM  images  of  thermally  induced  microdefects  in  Si 

54  MeV  electron  channeling  radiation  spectra  from  Si  (light  dots)  and  Si  +  O 

(heavy  dots) 

(100)  potentials  and  eigen  values  for  17  MeV  electrons  in  a)  GaAs,  b)Ala3 
Gap  7  As  and  c)  AlAs. 

Calculated  spectra  for  17  MeV  electrons  in  a)  GaAs,  b)  AIQ3  GaQ7As  and  c) 
AlAs. 


54  MeV  electron  channeling  radiation  spectra  from  Si,  a)  (100)  plane  at 
-190°C  (heavy  dots)  and  5°C  (light  dots),  b)  (110)  plane  at  -180°C  (heavy 
dots)  and  7°C  (light  dots). 

54  MeV  electron  channeling  radiation  peak  energy  vs.  temperature  in  Si.  a) 
(100)  and  b)  (110)  plane.  Solid  curves  show  theoretical  values. 

54  MeV  electron  channeling  radiation  linewidth  vs.  temperature  in  Si  a)  (100) 
and  b)  (1 10)  planes.  Dotted  lines  are  the  best  fit  curves. 

54  MeV  electron  channeling  radiation  spectra  from  (100)  planes  in  LiH  (light 
dots)  and  LiD  (heavy  dots) 

12.6  MeV  electron  channeling  radiation  spectra  from  Diamond  a)  <100>,  b) 
<1 10>,  and  c)  <1 1 1>  axis. 

Initially  simulated  isolation  devices  with  different  l-V  characteristics;  semi- 
recessed,  top,  and  fully-recessed,  bottom.  Tqx  =  6500  Angstroms,  =  0.5 

/Am. 


Cross-section  of  conventional  isolation  region  structure  with  a  plot  of  the 
surface  potential  versus  the  horizontal  position.  Tqx  =  6500  Angstroms, 
x.=  0.8 /tm. 

Cross-section  of  isolation  region  structure  with  the  corner  effect  and  a  plot 
of  the  surface  potential  showing  the  new  position  of  the  barriers  near  the 
oxide  corners.  TQx  =  6500  Angstroms,  =  0.25  /Am. 

Calculated  l-V  curves  for  the  devices  of  Figures  1-11  and  1-12.  The  device 
with  the  comer  effect  displays  better  isolation  than  the  normal  device 
without  the  corner  effect.  VDS  =  5  V,  vSUB  *  0  V. 

Plots  of  the  surface  potential  of  the  transistor  with  the  comer  effect  in  Figure 
1-12,  illustrating  the  formation  of  the  potential  barriers  at  the  corners  with 
increasing  gate  voltage.  Vs  =  0  V,  V0  =  0  V,  V$ug  =  0  V,  VQ  =  -0.8,  2,  5,  8,  11, 


14  V. 

PISCES  generated  plots  of  the  surface  potential  of  a  device  with  the  corner 
effect  showing  the  effect  of  drain  bias.  VQ  -  18  V,  Vs-  0  V,  VS(JB  ■■ 


=  ov,  v0  =  o, 


0.5, 3,  8  V.  Nsug  --2X  10,6cm'3. 

Cross-sections  of  isolation  region  structures  with  different  interface 
curvatures.  The  top  structure  has  a  larger  radius  than  the  bottom  structure. 

vos  =  5V- VSU8  =  0V- 

Calculated  I--V  curves  for  the  devices  of  Figure  1-16.  The  device  with  the 
sharp  corner  displays  better  isolation  than  the  device,  with  the  gradual 
corner. 

Cross-section  of  small  isolation  region  transistor  with  its  calculated  l-V 
characteristics.  The  gate  has  no  control  of  the  drain  current. 

Effects  of  substrate  bias.  The  cross-section  of  a  isolation  region  device  with 


8 

9 

12 

13 

15 

16 

17 

18 
20 
26 

27 

28 

29 

31 


32 

33 

34 

35 

36 


JSEP  Annual  Report 


Page  iii 


June  1984 


Stanford  University 


Stanford  Electronics  Laboratories 


Figure  1-20: 


Figure  1-21: 


Figure  1*22: 
Figure  1-23: 

Figure  1-24: 

Figure  1-25: 


Figure  1-26: 


Figure  1-27: 


Figure  1-28: 


Figure  1-29: 
Figure  1-30: 
Figure  1-31: 

Figure  1-32: 

Figure  1-33: 
Figure  1-34: 
Figure  1-35: 


the  corner  effect  shows  the  smaller  increase  in  the  depletion  width  at  the 
corner  than  at  the  center  planar  section.  A  transistor  with  the  corner  effect 
shows  a  larger  shift  in  the  calculated  l-V  curves  with  an  applied  substrate 
bias  than  one  without. 

Plot  of  the  subthreshold  slope  factor  versus  the  interface  radius.  The 
discrete  points  are  from  GEMINI  simulations,  with  the  one  on  the  far  right 
hand  side  from  standard  one  dimensional  analysis.  A  device  with  a  small 
interface  radius  will  have  a  flatter  subthreshold  slope  and  hence  better 
isolation. 

Cross  sections  of  a  test  transistor  during  the  fabrication  process.  Top:  after 
die  arsenic  has  been  implanted  in  the  source,  gate,  and  drain  regions. 
Upper  middle:  after  the  etching  of  the  silicon  in  the  gate  region.  Lower 
middle:  after  the  polysilicon  has  been  deposited  and  defined.  Bottom:  The 
completed  test  transistor. 

Microphotograph  of  experimental  test  structures. 

SEM  cross-section  of  isolation  structure  with  the  largest  fabricated  interface 
radius.  r=  1.00 pm. 

SEM  cross-section  of  isolation  structure  with  the  smallest  fabricated 
interface  radius,  r-  0.57  pm. 

Plot  of  measured  l-V  characteristics  of  four  devices  with  varying  amounts  of 
the  comer  effect.  Curve  A  is  from  a  device  with  no  corner  effect  Curves  B, 
C,  and  0  are  from  devices  with  decreasing  interface  radii.  VS(jg  •  0.0  V, 
Vds  =  0.05  V,  channel  length  is  7.5  pm,.  width  is  500  pm.  NSUBP  I  1.3  X 
70™  cm-3. 

Plot  of  measured  l-V  characteristics  of  two  devices  displaying  the  effect  of 
substrate  bias.  Curve  A:  no  corner  effect,  V$(J3  *  0.0  V;  Curve  B:  VS(JB  *  -0.5 
V.  Curve  C:  with  corner  effect,  VS(JB  *  0.0  V;  Curve  D:  VSUB  *  -0.5  V.  Channel 
length  is  2.5  pm,  width  is  500  pm.  There  is  a  larger  shift  in  the  curves  due  to 
applied  substrate  bias  in  the  device  with  the  comer  effect.  WCI1D«*  1.3  x 
1016  cm-3.  SUB 

Plot  of  measured  l-V  characteristics  of  two  devices  displaying  the  effect  of 
drain  voltage.  Curve  A:  no  comer  effect,  VDS  *  0.5  V;  Curve  B:  VQS  r  5.0 
V.  Curve  C:  with  corner  effect,  VQS  =  0.5  V;  Curve  D:  VQS  =  5.0  V.  Channel 
length  is  2.5  pm,  width  is  500  pm.  There  is  a  smaller  shift  in  the  curves  due 
to  the  applied  drain  voltage  in  the  device  with  the  corner  effect.  A/-..-  =  1.5 

4  OUtS 

X  1015  cm3. 

Plot  of  the  drain  current  versus  the  drain  voltage  with  a  constant  gate 
voltage.  Curve  A  is  from  a  device  with  no  comer  effect  and  VQ  =  -0.06 
V.  Curve  B  is  from  a  device  with  the  corner  effect  and  VQ  =  1.07  V.  Channel 
length  is  2.5  pm,  width  is  500  pm.  The  device  with  the  comer  effect  does  not 
experience  drain  induced  barrier  lowering.  Nsug=  1.5  x  I0,5cm'3. 
l-V  characteristic  of  an  MIS  device  with  64  angstrom  nitride 
l-V  characteristics  of  devices  with  40  angstrom  nitride 
Conduction  characteristics  of  devices  with  100  angstrom  oxide  and 
nitroxide 

High  frequency  C-V  for  MIS  devices  with  (a)  100  angstrom  oxide  and  (b)  60 
angstrom  thermal  nitride 

Breakdown  characteristics  of  devices  with  thermal  nitride  insulators 
Effect  of  trapping  on  F-N  conduction  characteristics 
Time  dependent  electrical  breakdown  characteristics  for  devices  with  nitride 
grown  at  1230  °C  for  (a)  1  hour  on  p-type  and  (b)  4  hours  on  n  type  Si 


37 


40 


41 

41 

42 

43 


44 


46 


47 


51 

52 

52 

53 

54 

55 

56 


JSEP  Annual  Report 


Page  iv 


June  1984 


Stanford  University 


Stanford  Electronics  Laboratories 


Figure  1  -36:  High-frequence  C-V  plots  for  different  finger  widths  (L),  or  perimeter  to  area  60 
ratiois  (P/A).  The  inversion  region  capacitance  is  proportional  to  P/A. 
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INTRODUCTION 


This  report  describes  work  performed  under  the  Joint  Services  Electronics  Program,  sponsored  by 
the  Departments  of  the  Air  Force  (Air  Force  Office  of  Scientific  Research),  Army  (Army  Research 
Office),  and  Navy  (Office  of  Naval  Research),  under  contract  DAAG29-81  -K-0057  at  the  Stanford 
Electronics  Laboratories,  Stanford  University,  Stanford,  California. 


The  three  most  significant  accomplishments,  as  determined  by  the  Director,  are  summarized  as: 
1.  Real-Time  Statistical  Data  Processing 


The  goals  of  this  research  program  include  focusing  on  methods  for  determining  the 
number  of  sources  when  die  signal  is  completely  coherent  with  one  or  more 
interferences.  A  new  approach  to  the  problem  of  determining  the  number  of  sources  for 
a  passive  array  was  developed  [7].  The  method  is  based  on  information  theoretic  criteria 
concepts  for  model  selection.  Rather  than  requiring  the  definition  of  thresholds  (as 
needed  in  the  conventional  hypothesis  testing  approaches  for  determining  the  number  of 
sources  in  the  eigenstructure  framework),  the  new  approach  does  not  require  any 
subjective  threshold  settings.  Simulation  studies  yield  good  results  under  varying 
situations  including  low  signal  to  noise  ratios,  small  number  of  data  samples  and  closely 
spaced  signals. 


The  new  approach  is  equally  applicable  in  other  areas  where  eigenstructure  methods 
have  been  used.  Some  examples  are  determining  the  number  sinusoids  in  time  series 
data,  the  number  of  poles  from  the  natural  response  of  a  linear  system  and  the  number  of 
overlapped  echoes  in  backscatter  data  typical  in  radar,  etc. 

2.  Physics  and  Technology  of  Submicron  Devices 

This  research  program  was  divided  into  a  number  of  projects.  Two  of  these  projects  have 
matured  to  the  point  that  these  seed  projects  have  spawned  new,  larger  projects  which 
endeavor  to  carry  these  studies  forward  and  to  study  the  application  of  the  capabilities 
shown  to  be  possible  with  the  work  sponsored  by  this  JSEP  contract. 


One  of  the  projects  has  been  aimed  at  developing  SIPOS  (Semi- Insulating  Poly  Silicon) 
emitters  for  bipolar  transistor  applications  [8, 6,  5].  Using  heterojunction  contacts  to 
emitter  regions  as  discussed  in  previous  JSEP  reports,  the  transistor  current  gain 
increased  by  an  order  of  magnitude  over  that  attainable  with  conventional  diffused 
emittters  and  metal  contacts.  This  makes  substantial  improvements  in  device 
performance  and  optimization  of  bipolar  devices  for  VLSI  applications  possible.  Follow- 
on  work  is  expect  to  continue  under  a  contract  under  negotiation  with  Sandia. 

In  another  of  the  projects,  a  mathematical  analysis  of  the  physics  of  trench-like  isolation 
region  structures  resulted  in  new  physical  understanding  of  these  structures  and  in  an 
equation  relating  the  leakage  current  to  the  radius  of  the  interface  [2, 1,  3].  The  strong 
dependence  of  the  corner  effect  on  the  radius  was  verified  by  the  fabrication  of  several 
test  transistors  with  different  radii.  Measurements  of  these  transistors  also  demonstrated 
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the  improved  isolation  properties  with  substrate  bias  and  drain  voltage.  Burying  the  oxide 
and  sharpening  oxide  corners  means  improved  isolation  and  more  compact  VLSI  layouts. 
This  work  is  continuing  under  DARPA  sponsorship  on  contracts  MDA903-79-C-0257  and 
MOA903-84-K-0062. 


3.  Signal  Processing  and  Compression 

This  research  program  was  also  divided  Into  a  number  of  projects.  Of  special  interest  is 
the  work  in  vector  quantization.  The  project  has  focused  on  three  principal  problems: 

The  design  of  feedback  vector  quantizers  for  data  compression,  the  design  of  adaptive 
vector  quantizers  for  speech  compression,  and  the  extension  of  vector  quantizer  design 
algorithms  to  the  design  of  combined  quantization/classification  distributed  sensor 
networks.  In  the  area  of  Feedback  Vector  Quantizers,  a  stochastic  gradient  algorithm 
for  designing  predictive  vector  quantizers  (PVQ)  was  developed  in  order  to  get  better 
codebooks  and  faster  design  speed.  To  further  improve  the  performance,  adaptive  VQ’s 
which  combine  the  waveform  coding  techniques  and  linear  predictive  coding  (LPC) 
techniques  were  studied. 

Codebook  and  predictor  design  based  on  a  stochastic  gradient  algorithm  simultaneously 
improves  the  linear  predictive  coefficients  and  the  vector  quantizer  in  the  feedback  loop 
iteratively  for  a  given  training  sequence.  This  provides  better  performance  and  faster 
convergence.  Further  details  are  reported  in  [4]. 

An  NSF  contract  will  be  supporting  the  related  theory  of  the  finite  state  vector  quantizers 
in  the  future.  Related  practical  work  is  now  being  pursued  at  the  University  of  California, 
Berkeley  and  at  Bell  Labs. 

This  report  contains: 

•  Summaries  of  work  for  die  JSEP-sponsored  projects  for  the  period  April  1,  1983  through 
March  31, 1984. 

•  A  List  of  JSEP  Sponsored  Publications  through  31  March  1984  is  included  at  the  end  of 
each  summary. 
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1.  VERY  LARGE  SCALE  INTEGRATION 

1.1  Properties  of  Materials:  Application  of  Channeling  Radiation 
to  a  Study  of  the  Properties  of  Materials 

Principal  Investigator:  R.  H.  Panted 

1.1.1  Scientific  Objectives 

Our  long  term  objective  is  to  determine  those  crystal  properties  which  can  be  studied  by  means  of 
channeling  radiation,  and  to  relate  these  properties  to  experimental  observations.  We  have 
successfully  applied  this  technique  to  a  variety  of  materials  (Si,  Ge,  LiF,  NaCI,  W,  and  LiH),  and  in  the 
future  the  emphasis  will  be  on  investigating  defects  and  properties  of  electronic  semiconductors,  (Si, 
GaAs,  and  AIGaAs).  These  materials  are  important  for  their  use  respectively  in  VLSI,  in  high  speed 
microwave  and  digital  components,  and  in  semiconductor  lasers. 

This  approach,  channeling  radiation  study  of  crystals,  is  nondestructive,  and  .  can  be  used  to 
investigate  electron  distribution,  crystal  potentials,  atomic  vibrations,  damage  and  defects,  and 
anistropic  properties.  In  comparison  to  the  electron  microscope  (TEM,  SEM  or  STEM),  channeling 
radiation  is  applicable  to  much  larger  and  thicker  samples  (up  to  100  jim  thick).  Channeling  radiation 
provides  an  alternative  technique  to  x-ray  diffraction  for  the  determination  of  crystal  parameters. 
Whereas  x-ray  measurements  rely  upon  the  accurate  measurement  of  diffracted  beam  intensities,  the 
channeling  radiation  technique  involves  a  measurement  of  spectral  peak  positions,  which  can  be 
determined  more  precisely. 

During  the  next  year  we  plan  to  perform  the  following: 

•  repeat  the  experiment  with  oxygen  platelets  in  silicon  and  investigate  the  statistics  of 
platelet  structures  and  their  properties.  Oxygen  in  Si  has  been  an  important  subject  for 
many  years  because  of  enhanced  gettering  of  metal  ions  in  Si  wafers  with  oxygen 
platelets. 

•  investigate  clustering  of  dopants  and  defects  in  both  AIGaAs  and  GaAs.  This  is  an 
interesting  subject  because  of  the  use  of  these  materials  in  fast  transistors  utilized  in 
advanced  computers  and  in  microwave  structures. 

•  study  the  effects  of  damage  induced  in  the  above  electronic  semiconductors  by  both 
particle  beams  and  electromagnetic  radiation  while  the  damage  is  occurring. 

•  investigate  the  above  electronic  semiconductor  materials  at  low  temperatures.  By 
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reducing  the  crystal  temperature,  we  can  study  the  various  characteristics  of  the 
materials  as  a  function  of  temperature,  such  as  thermal  vibration  amplitude,  electron 
scattering,  and  electron  distribution.  These  results  may  be  of  value  in  deciding  whether 
low  temperature  operation  of  VLSI  has  sufficient  advantages  to  justify  the  increased 
system  complexity. 

Channeling  radiation  provides  a  new  tool  for  investigating  the  properties  of  crystals.  We  have 
begun  to  utilize  this  tool  for  both  "perfect"  crystals  and  crystals  with  defects,  and  we  propose  to 
extend  our  work  to  new  materials  and  new  types  of  defects. 

1.1.2  Progress 

During  the  past  year  we  performed  measurements  of  channeling  radiation  from  54  MeV  electrons  in 
Si  with  an  oxygen  platelet  impurity;  calculated  potentials  and  eigenvalues  for  17  MeV  electrons  in 
GaAs,  AIGaAs,  and  AlAs;  measured  channeling  radiation  down  to  liquid  nitrogen  temperature; 
measured  and  compared  electron  channeling  radiation  form  LiH  and  LiD;  and  measured  12.6  MeV 
electron  channeling  radiation  from  Si,  diamond  and  LiF. 

1 .1 .2.1  Channeling  Radiation  from  Si  with  an  Oxygen  Platelet  Impurity 

Oxygen  precipitation  in  silicon  has  been  the  subject  of  many  studies  utilizing  different 
techniques.  [1 , 10]  Oxygen  is  the  main  impurity  present  in  Czochralski  (CZ)  grown  silicon  and  its 
presence  may  directly  or  indirectly  affect  the  electronic  properties  of  the  device. 

We  have  measured  54  MeV  electron  channeling  radiation  from  a  130  pm  thick  Si  crystal  with 
oxygen  platelets.  The  wafer  was  subjected  to  a  750°C  anneal  under  dry  02  for  100  hours,  and  the 
oxygen  concentration  in  this  crystal  was  measured  at  26  ppm.  This  crystal  shows  the  following 
defects  from  high  resolution  TEM  pictures. 

1.  Thin  (10-30  angstroms)  platelets  of  200-1500  angstroms  diameter  parallel  to  (100)  Si 
planes. 

2.  Small  stacking  faults  lying  in  (111)  planes  in  the  vicinity  of  the  (100)  planar  silica  defect.  A 
TEM  photograph  of  a  lattice  with  a  (100)  platelet  after  a  750°C  anneal  is  shown  in  Fig.  1-1 

The  54  MeV  electron  channeling  radiation  spectra  from  the  Si  +  0  crystal  are  shown  in  Fig.  1  -2,  and 
compared  to  those  from  a  40  pm  thick  Si  crystal  without  the  oxygen  platelets.  Figs.  l-2(a)-(c)  are  the 
spectra  from  (100),  (110),  and  (111)  planes,  respectively.  Heavy  dots  are  the  data  from  the  Si  +  0 
crystal,  and  light  dots  are  from  Si.  The  (100)  and  (111)  spectra  from  the  two  crystals  are  not  very 
different,  however,  the  (110)  plane  spectrum  from  the  Si  +  0  crystal  shows  large  linewidth  broadening 
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compared  with  that  from  Si.  Also,  there  are  peak  energy  down-shifts  in  1-*0  and  2-+1  transitions 
(Table  1-1).  The  transition  energy  changes  in  (110)  planes  are  probably  not  due  to  the  thickness 
differences  between  the  two  crystals,  because  this  effect  was  not  observable  in  (100)  and  (111) 
planes. 

We  previously  reported  the  downshift  of  radiation  energies  from  the  (100)  plane  in  diamond  with  a 
nitrogen  defect  [8].  This  is  a  reasonable  result  considering  the  distortion  in  the  (100)  plane  resulting 
from  the  platelet.  It  is  not  clear  why  we  observed  the  downshift  in  the  (110)  plane  in  the  Si  +  0  crystal 
rather  than  from  the  (100)  plane.  This  experiment  will  be  repeated  with  a  thinner  Si +  0  (30  micron) 
crystal  and  will  be  studied  further. 
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Plane 

(100) 

(110) 

(111) 


Table  1*1:  Comparison  of  Peak  Energy  and  Linewidth  from  Si  with  and  without 

Platelets 


Trans¬ 

Si 

(AO  pun) 

Si+0( 130um) 

ition 

Peak  Energy  Linewidth 

Peak  Energy 

Linewidth 

(keV) 

(keV) 

(keV) 

(keV) 

l-»0 

95.8±1.2 

17.4±0.6 

94.8±1.2 

14.6±0.4 

2— >1 

80.8±0.8 

10.3±0.4 

60.5±0.8 

10.6±Q.5 

3— *2 

37.7±1.2 

11.0±0.8 

37.8±1.0 

10.3±0.8 

l-»0 

121.7±l-6 

18.1±0.8 

116.1±2.0 

21.3±1.0 

2-+1 

87.1±1.2 

11.9±0.4 

86.4±1.2 

14.8±0.8 

3-+2 

62.8±1.0 

10.3±0.4 

82.1±1.2 

9.0±0.8 

4— »3 

48.0±0.8 

7.7±0.4 

48.8±1.0 

8.4;fc0.6 

5-44 

37.1±0.8 

7.1±0.4 

37.7±0.8 

9.0±0.8 

6-45 

28.6^0.8 

4.6±0.8 

28.8±0.8 

5.8±0.8 

Smal  1 
Peak 

95.0±1.0 

18.8±0.6 

93.6±1.0 

21.8±0.6 

Large 

Peak 

36.0±0.8 

25.6±0.3 

36.0±0.4 

29.0±0.3 
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1.1. 2.2  Calculated  Potentials  and  Eigenvalues  in  GaAs,  AIGaAs,  and  AlAs 

GaAs-AIGaAs  is  an  important  material  because  of  its  use  in  the  superfast  computers  [6]  and  in 
diode  lasers  [7]  for  fiber  communication.  There  have  been  some  interesting  studies  about  clustering 
of  Al  [4,  S]  and  strain  [2]  in  these  materials.  These  phenomena  should  be  easily  observable  by  both 
positron  and  electron  channeling  radiation. 

Figure  1-3  shows  the  calculated  planar  potentials  and  eigenvalues  for  17  MeV  electrons  channeled 
by  (100)  planes  in  GaAs,  Al03  GaQ7  As,  and  AlAs.  Fig.  1-4  shows  the  calculated  spectra,  assuming  20 
micron  of  crystal  thickness.  In  AIQ4  Gag  7  As,  if  the  Al  is  not  uniformly  distributed  or  is  clustered,  one 
will  not  observe  a  spectral  peak  and  perhaps  there  will  be  no  channeling  radiation  at  all. 


JSEP  Annual  Report 


Page  1 1 


June  1984 


JSEP  Annual  Report 


Stanford  University 


Stanford  Electronics  Laboratories 


1.1. 2.3  Low  Temperature  Channeling  Radiation 

We  measured  54  MeV  electron  channeling  radiation  from  a  Si  crystal  as  a  function  of  temperature. 
For  electron  channeling,  the  thermal  vibration  of  the  atoms,  which  causes  transitions  between  eigen 
states  as  well  as  dechanneling,  is  the  primary  contribution  to  tinewidth.  At  lower  temperature  the 
linewidth  is  reduced,  since  the  thermal  vibration  of  atoms  decreases  with  temperature.  Another  effect 
of  cooling  is  to  shift  the  spectral  lines  to  higher  frequencies.  The  potential  well  for  a  fixed  atom  is 
deeper  than  the  well  for  a  vibrating  atom,  and  in  ther  former  case  there  is  greater  separation  between 
energy  eigen  values. 

Figures  1-5(a)  and  (b)  show  the  54  MeV  electron  channeling  radiation  spectra  from  (100)  and  (110) 
planes  in  a  Si  crystal  at  room  temperature  and  at  a  low  temperature.  Heavy  dots  are  low  temperature 
data  (-190°C  for  (100)  plane  and  -160°C  for  (110)),  and  room  temperature  data  are  light  dots.  We  can 
observe  the  1— *0  transition  energy  shift  very  clearly.  Also,  some  line  narrowing  can  be  noticed. 
Figures  1-6(a)  and  (b)  show  the  measured  and  calculated  channeling  radiation  peak  energy  as  a 
function  of  temperature,  using  the  Debye  approximation  [3]  to  estimate  the  thermal  vibration 
amplitude. 

For  both  the  (100)  and  (110)  planes,  the  1-»0  transition  shows  the  greatest  discrepancy  between 
the  theoretical  calculation  and  the  experimental  result.  However,  the  higher  n-values  the  agreement 
is  better.  From  our  experimental  results,  we  can  deduce  a  Debye  temperature  of  475“K,  which  is 
somewhat  lower  than  the  value  given  in  the  literature  [543°K]  [3].  In  other  words,  our  experiment 
shows  larger  atomic  vibration  amplitudes  than  the  published  values  (0.085  angstroms  vs  0.075 
angstroms  at  room  temperature). 

Figures  i-7(a)  and  (b)  show  the  measured  linewidths  for  transitions  in  (100)  and  (110)  planes  as  a 
function  of  temperature.  The  dotted  lines  are  the  best  Fit  curves.  As  we  expected,  the  linewidth 
decreases  as  the  crystal  temperature  deceases.  The  effect  is  strongest  for  1-+0  transitions,  since 
these  states  lie  deepest  in  the  potential  well  and  have  the  most  overlap  with  the  vibrating  atoms. 

1.1. 2.4  Electron  Channeling  Radiation  from  LiH  and  LiD 

We  measured  54  MeV  electron  channeling  radiation  from  LiH  and  LiD  crystals.  This  is  the  first  time 
that  channeling  radiation  has  been  used  to  identify  isotopes  in  crystals.  Figure  1-8  shows  the 
measured  spectra  from  (100)  planes  in  LiH  and  LiD  crystals.  The  spectrum  from  LiH  is  in  light  dots 
and  heavy  dots  are  for  LiD.  Tbe  thicknesses  of  the  crystals  are  190  microns  and  350  microns  for  LiH 
and  LiD,  respectively.  From  the  experimental  results,  the  radiation  peak  energies  from  LiD  are  higher 


JSEP  Annual  Report 


Page  14 


June  1984 


Stanford  University 


Stanford  Electronics  Laboratories 


Figure  1  -5:  54  MeV  electron  channeling  radiation  spectra  from  Si, 

a)  (100)  plane  at -190°C  (heavy  dots)  and  5°C  (light  dots). 

b)  (110)  plane  at  -180° C  (heavy  dots)  and  7°C  (light  dots). 

than  LiH,  which  is  contrary  to  what  we  calculated.  Table  1-2  shows  the  calculated  and  measured 
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Figure  6  («) 


Figure  1  *6:  54  MeV  electron  channeling  radiation  peak  energy  vs.  temperature  in 
Si.  a)  (100)  and  b)  (110)  plane.  Solid  curves  show  theoretical  values. 

channeling  radiation  energies  from  LiH  and  LiO.  In  this  calculation,  the  published  thermal  vibration 
amplitudes  [9]  were  used. 

For  both  LiH  and  LiD,  the  calculations  do  not  fit  very  well;  the  calculated  values  are  too  high  for  LiH 
and  are  low  for  LiO.  The  most  significant  difference  is  that,  based  upon  theoretical  calculations,  the 
radiated  peak  energies  from  LiH  are  higher  than  those  from  LiD,  whereas  the  experimental  results 
show  the  opposite  trend. 

From  these  experimental  results,  we  can  conclude  that  the  atomic  thermal  vibrations  for  LiD  are 
smaller  than  those  for  LiH.  We  improved  the  fit  to  our  data  by  changing  the  thermal  vibration 
amplitudes  based  upon  minimizing  the  mean  square  error  between  theory  and  experiment.  Table 
1-3  gives  the  calculated  values  for  the  modified  thermal  amplitudes  along  with  the  measured  values. 
The  agreement  is  much  better  than  that  which  can  be  obtained  using  the  published  values  for 
vibrational  amplitudes. 
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Plgura  7  (b) 


■  Figu  re  1  •  7 :  54  MeV  electron  channeling  radiation  iinewidth  vs.  temperature  in 

Si  a)  (100)  and  b)  (110)  planes.  Dotted  lines  are  the  best  fit  curves. 

1.1. 2.5  12.6  MeV  Electron  Channeling  Radiation  from  Diamond,  Si  and  LIF 

■  This  experiment  was  done  to  test  whether  we  can  obtain  good  results  for  lower  energy  electrons 
from  the  accelerator  at  LLNL,  since  we  need  lower  beam  energies  to  analyze  high  Z  semiconductor 
materials  (e  g.  GaAs)  studies.  Figure  1-9  illustrates  the  axial  channeling  spectra  from  <100>,  <11 0> 
and  <11 1>  axes  in  diamond.  These  results  indicate  excellent  channeling  radiation  spectra,  and  prove 
that  the  beam  quality  is  comparable  or  better  than  that  for  the  30-60  MeV  electrons  we  have  used 
previously. 
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Figu  re  1  -S:  54  MeV  electron  channeling  radiation  spectra  from  (1 00)  planes  in 

LiH  (light  dots)  and  LiD  (heavy  dots) 
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Table  1-2: 

Comparison  of  the  Calculated  and  Measured  Energies  for  LiH  and  LiD 

1 

i- 

Calculated*(KeV) 

Measured  (KeV) 

r 

— = 

Plane 

Transition  L1H  LiD 

LiH 

LiD 

, 

1—0 

62.4  55.8 

57.03=0.6 

58.3=f0.4 

(100) 

2—1 

45.6  43.6 

41.83=0.4 

43.63=0.3 

■ 

3— >2 

33.1  32.9 

29.63=0.4 

30.33*0.3 

1 

4— *3 

23.6  23.7 

21.1=F0.5 

21.2=F0.4 

i 

(110) 

1-0 

47.7  42.7 

41.3=F0.6 

45.33=0.4 

2—1 

29.4  28.3 

26.3=F0.5 

28.13=0.3 

1-0 

53.6  47.0 

52.03=0.4 

52.83=0.4 

1 

(111) 

3—1 

32.6  31.5 

38.63=0.4 

37.4=j=0.3 

» 

2—1 

24.3  25.5 

28 . 0=p0 . 4 

25.13=0.3 

'Thermal  vibration  amplitudes  [9] 

R 

L1H  Li 

:  0.113  angstroms 

LID  Li: 

0.147  angstroms 

-  ; 

H:  0.153  angstroms 

D: 

0.178  angstroms 

Table  1-3: 

Comparison  of  Calculated  Energies  Using  Modified  Thermal  Amplitude 

’* 

and  Measured  Energies 

r 

i 

f 

Calculated* 

(Modified) 

Measured 

1 

Plane 

Transition  LiH  LID 

LiH 

LiD 

» 

1-0 

56.8  59.1 

57.03=0.6 

58.33=0.4 

(100) 

2—1 

42.3  44.2 

41.83=0.4 

43.63=0.3 

3—2 

31.5  32.8 

29.63=0.4 

30.33=0.3 

4—3 

23.0  23.7 

21.13=0.5 

21.23=0.4 

» 

(110) 

1-0 

43.1  45.1 

41.33=0.5 

45.33=0.4 

i 

2— *1 

26.7  28.3 

28.33=0.5 

28.13=0.3 

1-0 

51.6  52.4 

52.03=0.4 

52.83=0.4 

(111) 

3—1 

33.9  33.1 

36.5=F0.4 

37.43=0.3 

2—1 

30.0  28.5 

26.03=0.4 

25.13=0.3 

• 

•Thermal  vibration  amplitudes  used: 

L 

K 

LiH  LI: 

0.12  angstroms 

LiD  Li: 

0.12  angstroms 

I 

* 

H: 

0.24  angstroms 

D: 

0.20  angstroms 

1 

' 

- 

1 

1 

’) 

*1 
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Figure  1-9: 
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12.6  MeV  electron  channeling  radiation  spectra  from  Diamond  a) 
<100>,  b)<llO>,  and  c)  <1 1 1>  axis. 
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1 .2  PHYSICS  AND  TECHNOLOGY  OF  SUBMICRON  DEVICES 

Principal  Investigators:  C.R.  Heims,  R.F.W.  Pease,  J.D.  Plummer,  W.E.  Spicer,  R.M.  Swanson,  J.S. 
Harris 

Introduction 

This  research  program  is  divided  into  a  number  of  projects  which  are  described  separately  below, 
including  progress  in  each  area  during  the  past  year  and  an  outline  of  future  plans.  The  first  three, 
which  report  research  under  the  direction  of  J.D.  Plummer,  are  closely  connected  in  that  they  are 
aimed  at  submicron  VLSI  devices. 

In  any  MOS  or  bipolar  integrated  circuit,  there  are  three  generic  types  of  structures.  The  first  is  the 
active  switching/amplifying  device.  We  are  investigating  the  characteristics  of  alternative  gate 
dielectrics  for  submicron  NMOS  and  CMOS  circuits.  Thermally  grown  Si3N4  and  Si02/Si3N4 
composites  appear  to  have  significant  advantages  over  Si02.  This  first  project  is  aimed  at  evaluating 
material  properties  and  corresponding  device  characteristics  associated  with  this  alternative 
dielectric. 

The  second  generic  type  is  a  parasitic  device  and  is  the  result  of  the  local-oxidation  techniques 
commonly  used  for  device-to-device  isolation.  Such  isolation  is  generally  achieved  in  both  bipolar 
and  MOS  integrated  circuits  through  selective  field  oxidation  coupled  with  self-aligned  ion-implanted 
channel  stops.  Despite  its  current  technological  dominance,  this  isolation  technique  will  be 
inadequate  as  device  geometries  approach  1  pm  and,  ultimately,  submicron  sizes.  As  a  result,  the 
scientific  objectives  of  this  second  project  are  to: 

1.  use  two-dimensional  numerical  simulation  to  determine  and  better  understand  the  limits 
of  local-oxidation  device-isolation  schemes  and 

2.  investigate,  both  theoretically  and  experimentally,  isolation  techniques  that  may  improve 
the  isolation  characteristics  at  feature  sizes  of  <1  pm. 

The  third  generic  type  is  generally  a  passive  load  element  which,  in  bipolar  circuits,  often  takes  the 
form  of  diffused  or  ion-implanted  resistors.  In  NMOS  circuits,  the  most  common  device  is  a 
polysilicon  resistor.  Polysiiicon  layers  are  also  widely  used  as  gate  electrodes  and  for 
interconnections  in  MOS  circuits.  In  this  third  project,  we  are  investigating  the  limiting  behavior  of 
polycrystalline-silicon  layers  when  they  are  defined  with  lateral  and  vertical  dimensions  comparable  to 
the  grain  sizes  in  the  materian  The  specific  project  under  investigation  is  the  sensitivity  of  Si/Si02 
interface  properties  to  ambient  conditions  when  the  interface  is  located  under  a  polysilicon  electrode. 
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We  are  attempting  to  determine  the  extent  to  which  such  interfaces  are  sensitive  to  surrounding 
ambients  and  in  particular  have  looked  at  the  properties  of  interface  charges  in  small  20  structures. 

In  summary,  we  are  engaged  in  a  coherent  investigation  of  the  problems  of  submicron  silicon 
integrated  circuits,  and  the  three  projects  are  aimed  at  the  three  key  components  of  any  modem 
silicon  circuit.  Our  intent  is  to  establish  a  scientific  basis  for  understanding  and  circumventing 
problems  that  may  occur  in  scaling  Hie  three  generic  structures  to  submicron  dimensions. 

1.2.1  Electrical  Performance  and  Physics  of  Isolation  Region  Structures  tor  VLSI 

investigators:  S.  Goodwin,  J.  Plummer 

1.2. 1.1  Scientific  Objectives 

Virtually  all  modem  MOS  and  bipolar  integrated  circuits  use  the  local  oxidation  of  silicon  to  isolate 
laterally  between  devices  [1].  This  isolation  structure  is  typically  formed  by  doping  the  field  region 
using  ion  implantation  followed  by  the  growth  of  a  thick  oxide  in  the  field  region.  The  active  device 
regions  are  protected  from  these  steps  by  a  masking  layer  of  silicon  nitride  which  is  subsequently 
removed.  The  parasitic  device  that  is  formed  comprises  a  source  and  drain  from  the  source/drains  of 
two  adjacent  active  transistors  and  a  gate  made  of  an  interconnection  wire  of  polysilicon  or  metal. 
Good  isolation  entails  that  the  parasitic  device  be  turned  off  and  conduct  a  minimum  of  current, 
which  the  increased  doping  and  thicker  oxide  provide.  With  the  parasitic  transistor  conducting  a 
minimum  of  current,  there  is  less  crosstalk,  less  leakage  from  storage  nodes,  and  less  unwanted 
interaction  between  the  active  transistors.  This  technique  of  LOCOS  isolation  has  certain  advantages 
such  as:  the  self-alignment  of  the  field  oxide  to  the  more  heavily  doped  channel  stop,  ion  implantation 
of  the  channel  stop  to  allow  an  optimum  choice  of  substrate  doping  for  the  active  devices,  and  the 
reduction  of  step  height  between  the  field  and  active  device  regions. 

The  LOCOS  process  has  the  undesirable  by-product  of  a  "bird’s  beak”  due  to  lateral  oxidation, 
which  consumes  chip  area  and  therefore  reduces  device  packing  density.  Another  drawback  is  the 
lateral  diffusion  of  the  channel  stop  implant  during  the  field  oxidation  which  can  also  reduce  the 
packing  density  by  encroaching  on  active  device  channel  regions  and  thereby  requiring  larger  drawn 
device  dimensions.  Lateral  channel  stop  diffusion  also  degrades  the  source/drain  to  substrate 
capacitance.  Several  recent  papers  have  addressed  technological  alternatives  to  the  LOCOS 
process,  to  circumvent  these  problems.  The  proposed  alternative  technologies  can  be  divided 
roughly  into  two  groups.  One  group  is  a  return  to  the  previous  technique  of  a  planar 
oxide/semiconductor  interface  with  different  oxide  thicknesses  above  the  interface.  This  can  be 
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accomplished  by  growing  a  thick  field  oxide  everywhere  and  then  removing  the  oxide  in  the  active 
device  regions  [14],  or  by  depositing  polysilicon  and  selectively  oxidizing  it  into  the  field  oxide  [10]. 
Strict  attention  must  be  paid  to  controlling  the  step  up  to  the  field  regions  and  doping  of  the  field 
regions.  The  other  group  controls  the  lateral  growth  of  the  oxide  during  LOCOS  with  different 
processing  [4, 8, 9, 2, 13].  These  processing  variations  utilize  different  configurations  of  silicon 
nitride,  or  deposited  oxides.  We  have  investigated  in  this  work  some  of  the  implications  for  electrical 
performance  that  these  new  isolation  structures  imply.  By  analyzing  their  physics  of  operation,  the 
electrical  characteristics  can  be  predicted  for  new  generations  of  devices.  These  predictions  are 
experimentally  verified  through  the  fabrication  of  several  "trench-like"  structures.  Both  new 
technology  and  a  better  understanding  of  isolation  region  physics  will  be  necessary  to  solve  the 
isolation  problems  in  submicron  structures. 

1.2. 1.2  Progress 

Computer  Simulations 

To  obtain  a  preliminary  understanding  of  the  device  physics,  we  chose  to  use  computer  simulations 
for  flexibility  and  speed.  Recognizing  the  inherent  two-dimensional  nature,  the  abundance  of  non- 
planar  surfaces,  and  the  normal  operation  of  isolation  devices  in  the  subthreshold  regime,  we  used 
GEMINI,  a  program  that  solves  Poisson's  equation  in  two  dimensions  [7].  GEMINI  can  be  used  to 
generate  plots  of  the  potential  in  the  semiconductor  along  the  current  path.  This  aids  in  the 
understanding  of  the  device  operation  since  in  die  subthreshold  regime,  the  drain  current  is 
controlled  by  the  injection  of  carriers  over  a  potential  barrier  and  is  related  exponentially  to  die  height 
of  the  barrier. 

The  first  series  of  simulations  compared  two  variations  of  the  local  oxidation  technique,  a  semi- 
recessed  and  a  fully- recessed  structure  as  shown  in  Figure  1-10.  The  two  devices  were  identical, 
except  in  the  fully-recessed  structure  the  silicon  had  been  etched  in  the  field  regions  prior  to  the  field 
oxidation.  These  structures  are  somewhat  ideal  representations  of  actual  isolation  region  structures, 
and  were  chosen  to  simplify  the  initial  analysis.  More  realistic  structures  will  be  considered  later  in 
the  experimental  results  section  of  this  report.  The  two  structures  in  the  figure  have  the  same  oxide 
thickness,  the  same  doping  concentration,  and  the  same  junction  depths;  so  we  were  surprised  to 
find  distinctly  different  electrical  characteristics.  Comparing  their  calculated  log(/D)  vs.  VQS  curves, 
the  fully- recessed  device  conducted  less  current  because  of  a  flatter  subthreshold  slope,  and  hence 
had  better  isolation  properties. 

Examining  the  possible  causes  for  this  difference,  we  focussed  on  the  corner  of  the  field  oxide  and 
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Figu  re  1  •  1 0:  Initially  simulated  isolation  devices  with  different  I- V 

characteristics;  semi-recessed,  top,  and  fully- recessed,  bottom. 

Tox  ~  6500  AnQStroms,  =  0.5  /im. 

the  different  configurations  of  the  nearby  source  and  drain.  Extending  the  channel  of  an  isolation 
device  around  the  oxide  corner  provides  significant  improvement  in  isolation  due  to  geometrical 
reasons.  To  understand  this,  examine  a  normal  device  whose  source  and  drain  junctions  extend 
below  the  oxide/semiconductor  interface,  as  shown  in  Figure  1-11.  Here  the  channel  does  not 
include  the  comers  and  standard  one-dimensional  theory  applies.  As  the  calculated  plot  of  surface 
potential  along  the  Si/Si02  interface  shows,  there  is  a  potential  barrier  centered  under  the  gate.  As 
the  gate  voltage  is  increased,  the  height  of  the  barrier  decreases,  and  more  current  is  able  to 
surmount  the  barrier  and  reach  the  drain. 

Figure  1-12  shows  the  identical  device  except  the  source  and  drain  junctions  are  now  much 
shallower,  lying  above  the  oxide/semiconductor  interface.  By  pulling  these  junctions  up,  the  oxide 
corners  have  been  "exposed"  and  incorporated  into  the  channel.  Now  the  current  flow  and  device 
physics  become  very  two-dimensional  in  nature.  One-dimensional  theory  fails  to  hold,  as  seen  by  the 
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VG»10  Vs=0  V0»0  VSUB  =  0  NA  =  2  x  1016 cm'3 

Figu  re  1  •  1 1 :  Cross-section  of  conventional  isolation  region  structure 
with  a  plot  of  the  surface  potential  versus  the  horizontal 
position.  Tox  -  6500  Angstroms,  xf  -  0,8  /*m. 

calculated  surface  potential  of  the  semiconductor,  in  the  comer  regions  of  the  thick  field  oxide. 
Instead  of  the  barrier  being  centered  under  the  gate,  it  is  now  located  at  the  oxide  corners.  This  can 
be  understood  by  examining  the  two  crosshatched  regions  which  are  bounded  by  the 
oxide/semiconductor  interface,  the  depletion  edge,  and  two  electric  field  lines.  At  the  planar  region 
in  the  middle  of  the  channel,  this  region  is  rectangular  in  shape,  while  the  analogous  region  at  the 
comer  of  the  oxide  is  pie-shaped.  Thus  for  equal  interface  areas  and  depletion  widths,  the  pie¬ 
shaped  region  will  contain  more  bound  depletion  charge  than  is  needed  to  support  the  gate  voltage; 
therefore  the  depletion  width  is  less  at  the  corner  than  at  the  middle  of  the  channel.  The  depletion 
width  is  thinner  because  the  depletion  charge  supported  by  the  gate  voltage  is  more  spread  out. 
Because  the  depletion  width  is  smaller,  there  is  less  band  bending  and  hence  the  potential  at  the 
comer  is  less  as  well.  This  results  in  the  surface  potential  distribution  shown  in  Figure  1-12.  As  the 
gate  voltage  increases  from  a  low  value,  the  potential  at  the  comers  increases  much  more  slowly  than 
in  the  middle  of  the  channel. 
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VG  *  18  Vs=0  Vq  =  0  VSUB  =  0  NA=  2  x  I0,6cm'3 

Figu  re  1  •  1 2 :  Cross-section  of  isolation  region  structure  with  the  corner 

effect  and  a  plot  of  the  surface  potential  showing  the  new  position  of 
the  barriers  near  the  oxide  comers.  Tox  =  6500  Angstroms, 

Xj  =  0.25  pm. 

The  two  devices  of  Figures  1-11  and  1-12  have  been  biased  in  the  simulations  to  produce  the  same 
drain  current  of  4pA//*m  of  channel  width.  Because  of  the  corner  effect,  the  gate  voltage  is  10  V  in 
Figure  1-11  but  is  increased  to  18  V  for  the  device  in  Figure  1-12.  The  higher  gate  voltage  has 
produced  a  strong  inversion  condition  in  the  center  region  of  Figure  1-12;  the  current  is  then  limited 
by  the  potential  barriers  at  the  corner  regions. 

A  second  way  of  understanding  the  corner  effect  is  to  recognize  that  at  the  corner,  the  electric  field 
lines  diffuse  and  spread  out  as  they  penetrate  into  the  semiconductor.  Thus  there  is  a  weaker  link  of 
the  potential  in  the  semiconductor  to  the  gate  voltage,  as  shown  by  the  plot  of  the  potential.  This 
viewpoint  reveals  that  the  corner  effect  is  not  limited  to  the  semiconductor  but  also  occurs  in  the 
oxide  because  of  its  shape.  The  electric  field  lines  emanating  from  the  gate  begin  spreading  in  the 
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oxide  and  continue  to  do  so  in  the  semiconductor.  The  spreading  or  fringing  of  the  electric  field  lines 
and  its  effect  on  electrical  characteristics,  is  similar  to  the  narrow  width  effect  seen  in  the  shrinking  of 
active  devices. 

The  different  characteristics  of  the  corner  device,  which  is  essentially  a  modified  concave  MOSFET 
[11],  allows  the  modeling  of  an  isolation  device  by  three  devices  whose  channels  are  in  series,  two 
corner  devices  and  a  planar  device.  Because  of  the  geometrical  effects  explained  above,  the  corner 
device  has  a  flatter  subthreshoid  slope  and  a  higher  threshold  voltage,  and  thereby  controls  the 
current  flowing  in  the  isolation  device.  Figure  1-13  is  a  plot  of  the  calculated  log(/Q)  vs.  V'QS  curves  for 
the  devices  depicted  in  Figures  1-11  and  1-12.  The  device  with  the  corner  effect  clearly  has  better 
isolation  due  to  the  flatter  subthreshold  slope  as  well  as  the  higher  threshold  voltage.  (Note  that  the 
most  useful  definition  of  "threshold  voltage”  in  these  devices  is  probably  the  gate  voltage  necessary 
to  induce  a  leakage  or  subthreshold  current  of  some  value  —  perhaps  1  pA  per  micron  of  device 
width.  Using  this  definition  of  threshold  voltage,  the  corner  effect  increases  it  from  9  V  to  17  V.)  Thus 
for  a  given  gate  voltage,  the  device  with  the  corner  effect  passes  much  less  current  than  the  device 
without. 


A.  NORMAL  DEVICE 
B.  DEVICE  WITH  CORNER  EFFECT 

Figure  1-13:  Calculated  I- V  curves  for  the  devices  of  Figures  1-11  and  1-12. 

The  device  with  the  corner  effect  displays  better  isolation  than  the  normal 
device  without  the  corner  effect. 

vos  =  5V-W0V- 


The  corner  effect  requires  an  applied  gate  voltage  to  exist  and  affect  the  device  physics.  Initially 
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when  the  potential  in  the  semiconductor  is  unperturbed  by  the  gate  voltage,  there  is  no  spatial 
variation  in  the  surface  potential  between  the  comer  and  central  planar  regions.  Only  when  the  gate 
voltage  is  increased  and  the  semiconductor  surface  is  perturbed  by  it,  does  the  corner  effect  appear. 
Figure  1-14  shows  several  calculated  plots  of  the  surface  potential  along  the  Si/Si02  interface  for 
several  different  voltages  applied  to  the  gate  of  the  device  shown  in  Figure  1-12.  Notice  that  the 
comer  effect  does  not  emerge  until  the  gate  voltage  begins  to  perturb  the  surface  potential.  As  a 
result,  a  transistor  with  the  comer  effect  goes  through  several  stages  of  operation  as  the  gate  voltage 
is  increased.  At  low  gate  voltages,  the  surface  potential  is  only  slightly  affected  by  the  gate  and  the 
transistor  behaves  in  a  conventional  manner  with  the  entire  channel  acting  as  a  barrier  to  current 
injection  from  the  source.  The  drain  current  then  shows  a  l/L  dependence  because  the  entire 
channel  contributes  to  the  potential  barrier.  As  the  gate  voltage  is  increased  the  transistor  operation 
enters  a  transition  zone  where  the  potential  barrier  is  slightly  larger  at  the  oxide  corners  than  at  the 
central  planar  region.  The  planar  region  thus  has  a  diminished  effect  on  the  current  flow  and 
transistors  operated  in  this  regime  have  a  weak  1/L  dependence.  When  the  gate  voltage  is  further 
increased,  the  transistor  operates  in  a  mode  where  the  current  flow  is  controlled  only  by  the  comers 
and  is  unaffected  by  the  center  region.  This  is  because  the  potential  barrier  is  much  larger  at  the 
corners  than  in  the  middle  of  die  channel.  When  the  comer  device  dominates  the  parasitic  transistor, 
the  effective  channel  length  of  the  parasitic  transistor  changes  from  the  distance  between  the  source 
and  drain,  to  the  distance  along  the  curved  portion  of  the  oxide/semiconductor  interface.  As  a  result, 
transistors  with  the  comer  effect  and  different  source  to  drain  spacings,  do  not  show  the  typical  l/L 
dependence  in  subthreshold  current,  but  are  approximately  independent  of  L. 

For  the  potential  barrier  to  move  from  die  center  of  the  channel  to  the  corner,  the  channel  must 
exist  around  the  comer  of  the  oxide,  and  be  controlled  by  the  gate.  This  means  not  only  the 
metallurgical  boundary  of  the  source,  but  also  the  portion  of  the  depletion  region  controlled  by  the 
source,  must  lie  above  the  corner  for  the  isolation  to  be  improved.  Thus  the  junction  depth  is  of 
crucial  importance  to  the  comer  effect.  The  position  of  the  potential  barrier  at  the  comer  implies  that 
the  device  can  be  shortened  without  degrading  the  isolation  properties  of  the  device.  The  cause  of 
short  channel  effects,  drain  induced  barrier  lowering,  is  relatively  ineffective  due  to  the  new  location 
of  the  potential  barrier.  The  drain  voltage  must  affect  the  surface  potential  at  the  oxide  corner  near 
the  source  to  alter  the  controlling  potential  barrier.  If  the  electric  field  lines  from  the  drain  do  reach 
the  oxide  corner,  their  effect  is  further  reduced  by  the  same  geometrical  arguments  that  limit  the 
influence  of  the  gate. 

Another  computer  program  was  used  to  analyze  the  effect  of  drain  biases  on  the  device  physics  of 
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Figu  re  1  •  1 4:  Plots  of  the  surface  potential  of  the  transistor  with  the 
corner  effect  in  Figure  1-12,  illustrating  the  formation  of  the 
potential  barriers  at  the  corners  with  increasing  gate  voltage. 

Vs*0V,Vo*0V,  Vsua  *  0  V,  VQ  r  -0.8.  2,  5,  8.  11.  14  V. 

the  isolation  transistor.  The  GEMINI  program  has  limited  usefulness  for  understanding  structures 
with  the  comer  effect  and  applied  drain  biases  due  to  its  treatment  of  inversion  layers.  The  computer 
program  PISCES  solves  Poisson’s  equation  and  the  one  carrier  continuity  equation  in  two 
dimensions,  for  structures  with  non-planar  interfaces  [12].  This  enables  PISCES  to  simulate  parasitic 
transistors  with  the  corner  effect,  in  the  subthreshold  regime  with  appreciable  drain  biases,  and  also 
in  the  regime  where  the  device  is  turned  on.  Figure  1-15  contains  several  surface  potential  plots  of 
the  device  illustrated  in  Figure  1-12  for  different  drain  biases.  These  plots  illustrate  that  the  potential 
barrier  adjacent  to  the  source  is  not  affected  by  the  drain  potential. 

The  potential  barrier  adjacent  to  the  drain,  however,  is  modulated  by  the  drain  potential.  At 
sufficiently  high  drain  biases,  the  electric  field  from  the  drain  supports  a  large  amount  of  bound 
depletion  charge  around  the  adjacent  oxide  corner.  This  results  in  more  band  bending  there,  a 
reduction  in  the  height  of  the  potential  barrier,  and  eventual  elimination  of  the  barrier  adjacent  to  the 
drain.  The  removal  of  the  drain  potential  barrier  leaves  only  the  potential  barrier  adjacent  to  the 
source  providing  the  isolation.  If  the  gate  voltage  has  inverted  the  semiconductor  in  the  middle  of  the 
channel,  the  reduction  in  the  height  of  the  drain  barrier  allows  the  pool  of  minority  carriers  to  escape, 
which  changes  the  channel  into  a  state  of  deep  depletion. 
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Figure  1-15:  PISCES  generated  plots  of  the  surface  potential  of  a  device  with 
the  corner  effect  showing  the  effect  of  drain  bias.  vQ  =  18  V, 
Vs°0\r,VSUB*0V,V  :0,05.3.8V.  Nsub  =  2X  . 

1016  cm"3. 

The  curvature  of  the  oxide/semiconductor  interface  plays  a  major  role  in  the  corner  effect.  By 
varying  the  radius  of  the  interface  as  in  the  devices  of  Figure  1-16,  very  large  differences  in  electrical 
characteristics  result,  as  seen  in  the  plots  of  the  calculated  subthreshold  currents  in  Figure  1-17  for 
those  two  devices.  The  gate  voltage  necessary  for  IpA/pm  of  channel  width  increases  from  11  V  to 
37  V.  A  normal  planar  device  is  the  limit  of  a  corner  device  with  a  very  large  radius.  A  device  with  a 
radius  of  zero,  or  a  sharp  corner,  has  a  very  large  corner  effect  due  to  the  difficulty  the  electric  field 
experiences  in  reaching  into  the  corner  and  the  rapid  increase  in  available  depletion  charge  as  the 
depletion  width  increases.  Thus  corner  effect  devices  with  sharper  corners  provide  better  isolation  by 
having  flatter  subthreshold  slopes  and  higher  threshold  voltages.  The  angle  of  the  oxide  corner  also 
plays  a  role  in  determining  the  amount  of  the  corner  effect.  The  lower  device  of  Figure  1-16  has  a 
radius  of  zero  and  a  right  angle  corner.  As  the  angle  increases,  the  electric  field  will  be  able  to 
penetrate  the  corner  better  and  the  corner  effect  will  be  lessened.  When  the  angle  approaches  180° 
the  corner  effect  will  disappear. 

A  possible  very  small  isolation  structure  resulting  from  the  application  of  these  principles  is  shown 
in  Figure  1-18.  The  channel  ‘length’  has  been  shortened  to  0.2  pm  and  the  oxide  thickness  increased 
to  0.95  jim.  The  additional  solid  lines  in  the  figure  are  equipotential  curves.  The  lower  portion  of 
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Figure  1-16:  Cross-sections  of  isolation  region  structures  with 

different  interface  curvatures.  The  top  structure  has  a  larger  radius 
than  the  bottom  structure.  VQS  =  5  V,  VS{/8  =  0  V. 

Figure  1-18  shows  a  plot  of  the  calculated  log{/0)  vs.  VQS  curve  for  the  device.  This  transistor  has  a 
flat  subthreshold  current  curve  because  the  gate  has  no  control  of  the  surface  potential  at  the  comer, 
which  can  be  seen  by  the  shape  of  the  calculated  equipotential  curves.  The  drain  controls  the 
potential  barrier.  The  drain  current  that  flows  is  a  function  of  the  voltage  on  the  drain,  the  distance 
from  the  drain  to  the  corner,  and  the  amount  of  the  corner  effect.  Thus  the  current  can  be  controlled 
by  the  depth  of  the  oxide  cut;  the  deeper  the  cut  —  the  less  current.  This  is  because  the  depth  of  the 
cut  determines  which  equipotential  curve  intersects  the  lower  left  corner  of  the  oxide,  and  hence  the 
potential  barrier  height.  A  structure  such  as  this  could  have  considerable  use  in  the  isolation  between 
adjacent  diffused  lines  in  a  chip,  by  allowing  a  much  closer  spacing  of  them. 

The  addition  of  a  substrate  bias  enhances  the  corner  effect.  The  additional  gate  voltage  required 
for  the  extra  band  bending  changes  the  surface  potential  less  at  the  corner  than  in  the  planar  region 
and  hence  heightens  the  corner  effect.  Figure  1-19  shows  the  smaller  increase  of  depletion  width  at 
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A.  GRADUAL  CORNER  B.  SHARP  CORNER 

Figure  1-17:  Calculated  I--V  curves  for  the  devices  of  Figure  1-16.  The 
device  with  the  sharp  comer  displays  better  isolation 
than  the  device  with  the  gradual  comer. 

the  comer  for  an  applied  substrate  bias  of  -2  V  and  the  greater  shift  obtained  in  the  calculated  Iog(/D) 
vs.  VGS  curves  for  devices  with  the  comer  effect.  The  subthreshold  curves  of  the  device  with  the 
corner  effect  (curves  B  and  D)  shift  25  V  whife  the  curves  of  the  device  without  (A  and  C)  shift  13  V. 

Changes  in  the  doping  of  the  field  transistor's  channel  affect  the  corner  effect  in  different  ways,  so 
that  the  overall  change  in  electrical  characteristics  will  depend  on  the  configuration  of  the  particular 
structure.  Heavier  doping  results  in  the  depletion  edges  of  the  source  and  drain  being  pulled  up  and 
thereby  producing  a  greater  corner  effect  by  including  more  of  the  corner  in  the  channel.  Heavier 
doping  also  results  in  a  thinner  depletion  region  around  the  corner  which  reduces  the  two- 
dimensional  aspect  of  the  corner  in  the  semiconductor  and  lessens  the  comer  effect.  From  a 
practical  point  of  view,  only  a  small  range  of  substrate  doping  is  of  interest  for  isolation  structures, 
since  to  minimize  junction  capacitances,  the  smallest  doping  consistent  with  the  required  threshold 
voltage  would  be  used.  The  transistors  in  this  report  all  have  uniformly  doped  substrates,  but  the 
comer  effect  also  applies  to  transistors  with  localized  field  implants.  If  the  device  is  fabricated  in  such 
a  way  that  the  field  implant  extends  out  to  the  comers,  the  corner  effect  will  still  dominate  the  device 
physics.  Simulations  of  isolation  devices  with  field  implants  on  lightly  doped  substrates  have  shown 
that  dopant  concentration  can  be  less  at  the  corners  than  in  the  middle  of  the  transistor,  and  yet  the 
corner  effect  can  still  enhance  the  isolation.  This  is  because  the  geometrical  weakening  of  the 
linkage  between  the  gate  and  the  substrate  outweighs  the  lighter  doping  in  the  corner. 
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VGS  (volts) 

Figu  re  1  - 1 8:  Cross-section  of  small  isolation  region  transistor  with  its 

calculated  l-V  characteristics.  The  gate  has  no  control  of  the  drain 

current. 

Analytical  Development 

The  previous  section  has  demonstrated  with  the  aid  of  specific  computer  simulations,  the  basic 
ideas  which  can  be  used  in  optimizing  isolation  region  structures.  When  actually  optimizing  these 
structures,  however,  it  would  be  a  great  value  to  have  a  reasonably  simple  analytic  expression  which 
relates  the  electrical  properties  of  the  isolation  structure  (threshold  voltage,  subthreshold  slope,  etc.) 
to  the  physical  properties  of  the  structure  (trench  depth,  oxide  thickness,  corner  radius,  etc.).  Such 
an  analytical  relationship  would  be  faster  and  more  powerful  than  the  computer  simulations  and 
would  also  aid  in  estimating  the  sensitivity  of  the  electrical  properties  to  structural  changes  as  well. 

A  first  order  analysis  of  the  corner  effect  can  be  accomplished  by  considering  the  physics  of  the 
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A.  NO  CORNER  EFFECT  VSuB*0 
8.  CORNER  EFFECT  Vsua  *  0 
C.  NO  CORNER  EFFECT  VSUB  *  -2 
0.  CORNER  EFFECT  VSUB  =  -2 


Figure  1-19:  Effects  of  substrate  bias.  The  cross-section  of  a 

isolation  region  device  with  the  corner  effect  shows  the  smaller 
increase  in  the  depletion  width  at  the  comer  than  at  the  center 
planar  section.  A  transistor  with  the  corner  effect  shows  a  larger 
shift  in  the  calculated  I--V  curves  with  an  applied  substrate 
bias  than  one  without 


corner  device  and  assuming  the  corner  device  dominates  the  parasitic  transistor.  This  has  been  done 
in  two  separate  parts  [6].  The  first  part  determined  the  capacitance  of  the  oxide  structure  with 
conformal  mapping.  This  produces  the  variation  of  the  capacitance  along  the  oxide/semiconductor 
interface  in  the  comer  region  with  the  geometry  of  the  oxide.  The  second  portion  was  the  derivation 
of  the  MOS  equations  in  the  semiconductor  comer  region.  This  was  done  by  solving  the  one¬ 
dimensional  equations  in  cylindrical  coordinates.  Then  the  two  parts  were  brought  together  for  the 
complete  analysis  of  the  corner  effect. 


The  details  of  the  analytic  development  are  contained  in  [6].  We  show  here  only  an  example  of  the 
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results.  Figure  1-20  shows  a  plot  of  the  parasitic  transistor  inverse  subthreshold  slope  n  versus  r  for 
nsub  ' 2  x  1q16  cm'3,  T  -  6500  Angstroms,  and  <ps  =  0.5V.  The  plot  shows  a  dramatic  increase  in  n  as 
the  radius  decreases.  This  translates  into  greatly  improved  isolation  because  the  semiconductor  is 
less  influenced  by  the  gate  for  very  small  r.  Also  shown  in  Figure  1-20  are  several  discrete  points  that 
correspond  to  subthreshold  slope  data  acquired  from  GEMINI  simulations.  The  point  at  the  far  right 
hand  side  of  the  plot  is  from  the  one-dimensional  analysis,  or  a  radius  of  infinity.  The  point  for  r  =  o 
has  a  finite  slope  factor  instead  of  the  infinite  value  predicted  by  the  analytical  development.  This  is 
due  to  simplifying  assumptions  in  the  analytic  development. 


Figure  1  *20:  Plot  of  the  subthreshold  slope  factor  versus  the  interface 

radius.  The  discrete  points  are  from  GEMINI  simulations,  with  the  one 
on  the  far  right  hand  side  from  standard  one  dimensional  analysis.  A  device 
with  a  small  interface  radius  will  have  a  flatter  subthreshold  slope 
and  hence  better  isolation. 
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Figure  1-20  reveals  a  direct  link  between  an  electrical  parameter,  the  subthreshold  slope,  and  a 
structural  parameter,  the  corner  radius,  it  is  clear  that  subthreshold  slope  is  sensitive  to  the  value  of 
the  comer  radius  and  becomes  very  sensitive  for  values  below  0.25  pm.  Thus  in  designing  a  process 
for  the  fabrication  of  isolation  structures  for  VLSI,  one  would  want  a  small  corner  radius  in  order  to 
provide  good  isolation.  This  influences  the  choice  of  how  a  trench  is  etched  in  the  silicon  substrate 
and  also  how  it  is  filled.  If  the  process  depends  on  the  sharp  comers  to  provide  the  necessary 
isolation,  great  care  must  be  taken  to  ensure  that  the  variations  in  the  radius  are  sufficiently  small  to 
minimize  leakage  currents  in  worst  case  processing  variations.  The  analysis  can  be  expanded  to 
provide  an  equation  for  the  subthreshold  current  as  a  function  of  the  gate  voltage  and  structural 
parameters  [6]. 

Smaller  radii  and  sharper  angles  improve  the  isolation  according  to  the  device  physics  of  the 
structure’s  geometry.  That  geometry,  however,  may  lead  to  an  increased  number  of  crystal  defects 
because  of  stresses  due  to  the  fabrication  [5].  These  defects  can  cause  excessive  leakage  currents 
and  hence  alter  the  electrical  characteristics  of  the  device.  A  fabrication  process  must  consider  not 
only  the  geometry  of  the  isolation  device,  but  also  the  material  qualities  that  result  from  that  structure, 
in  arriving  at  an  optimal  structure. 

Experimental  Results 

Devices  have  been  fabricated  to  verify  the  comer  effect  experimentally.  Because  of  the  very  low 
current  levels  of  the  subthreshold  regime,  the  devices  were  designed  with  large  channel  widths  of  500 
pm  and  2820  pm.  Channel  lengths  varied  from  1 .25  to  40  pm  on  the  mask.  In  initial  fabrication  runs, 
the  devices  were  made  on  <1 10>  oriented  silicon  wafers  with  a  KOH  etch  to  produce  vertical  side  wall 
grooves  that  were  later  oxidized  to  form  the  field  oxide  isolation  devices.  That  orientation  of  silicon 
was  chosen  to  obtain  sharp  right  angle  comers  in  the  grooves.  The  bottom  surfaces  of  the  etched 
grooves  were  not  as  smooth  as  desired  however,  and  to  avoid  such  problems,  the  devices  described 
here  were  made  on  <100>  wafers  with  plasma  etching  of  the  silicon  to  produce  similar  structures. 
Arsenic  was  used  for  the  source  and  drain  to  ensure  shallow  junctions  above  the 
oxide/semiconductor  interface.  The  relative  position  of  the  junction  depth  and  oxide  interface  was 
varied  by  changing  the  junction  drive-in  time,  to  produce  structures  with  and  without  corner  effects. 
Different  oxide/semiconductor  interface  radii,  r,  were  produced  by  controlling  the  silicon  etching  and 
subsequent  oxidation. 

The  transistors  were  fabricated  on  p-type  6  0-cm  wafers.  A  portion  of  each  wafer  was  implanted 
with  boron  (100  keV,  9  x  1012  cm'2)  to  produce  other  transistors  with  a  heavier  substrate  doping. 
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Subsequent  high  temperature  processing  spread  this  implant  out  enough  to  consider  the  doping  to 
be  uniform.  Figure  1-21  displays  some  cross  sections  of  a  test  transistor  during  the  processing 
sequence.  Arsenic  was  implanted  (120  keV,  7  x  1015  cm'2)  in  the  source,  drain,  and  gate  regions.  A 
pad  oxide  was  then  grown  and  silicon  nitride  deposited.  Following  patterning  using  photoresist,  the 
exposed  Si3N4  and  Si02  were  removed  from  the  regions  in  which  the  isolation  regions  were  to  be 
formed.  This  was  followed  by  the  etching  of  the  silicon  substrate,  still  using  the  same  patterned 
photoresist  as  a  mask.  By  etching  through  the  arsenic  layer,  the  source  and  drain  were  defined.  The 
gate  region  was  locally  oxidized  using  the  Si3N4  as  a  mask  over  the  source  and  drain.  The  remainder 
of  the  processing  was  a  standard  poiysilicon  gate  NMOS  process.  The  resulting  test  structures  are 
shown  in  Figure  1-22. 

Three  different  interface  radii  were  fabricated  to  investigate  varying  amounts  of  the  corner  effect 
The  structure  with  the  largest  oxide/semiconductor  interface  radius  was  fabricated  by  isotropically 
etching  a  trench  and  then  thermally  growing  a  6400  Angstroms  oxide  at  950°C  in  H20.  Plasma 
etching  with  SFfi  was  used  to  produce  curved  sidewalls  whose  radii  were  increased  further  by  the 
oxidation.  Figure  1-23  contains  an  SEM  cross-section  of  this  structure.  An  intermediate  radius 
structure  was  fabricated  by  anisotropically  etching  the  trench  and  then  growing  a  6400  Angstroms 
oxide.  The  use  of  Freon  115  (C2CIFg)  for  the  plasma  etching  produced  vertical  sidewalls  with 
relatively  sharp  comers.  The  subsequent  oxidation  substantially  increased  the  radii  of  the  oxide 
comers  because  of  the  20  oxidation  process  in  the  comer  regions.  Anisotropic  etching  of  the  trench 
with  Freon  115,  followed  by  the  growth  of  3000  Angstroms  of  oxide  and  then  the  deposition  of  3400 
Angstroms  more  resulted  in  the  third  device  structure  which  had  the  smallest  interface  radius.  An 
additional  inert  ambient  anneal  was  done  for  these  structures  to  compensate  for  the  shorter  oxidation 
time.  Figure  1-24  contains  an  SEM  cross-section  of  this  structure.  The  radii  of  the  corners  produced 
by  the  Freon  115  did  not  increase  as  much  because  of  the  thinner  thermal  oxide.  The  radius  of  the 
interface  is  substantially  smaller  in  Figure  1-24  than  in  Figure  1-23,  0.57  jim  compared  to  1.00  /im. 
The  comer  radius  of  Figure  1-24,  however,  is  significantly  larger  than  the  radius  produced  by  the 
Freon  115  etching  because  of  the  2D  effects  of  the  thermal  oxidation.  It  should  also  be  noted  that  in 
all  of  these  structures  the  polysilicon  gate  lies  below  the  original  silicon  surface  since  the  etched 
trenches  were  not  completely  refilled  by  the  oxidation  and/or  Si02  deposition  processes.  More  ideal 
isolation  region  structures  (similar  to  Figure  1-12)  could  be  obtained  with  a  more  advanced 
technology  which  provided  complete  refilling  and  planarization. 

Figure  1  -25  shows  the  subth/eshold  characteristics  of  four  typical  devices.  All  have  channel  widths 
of  500  pm  and  channel  lengths  of  7.5  jim  on  the  mask.  The  oxide  thickness  is  6400  Angstroms  and 
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Figu  re  1  - 2 1 :  Cross  sections  of  a  test  transistor  during  the  fabrication 

process.  Top:  after  the  arsenic  has  been  implanted  in  the  source, 
gate,  and  drain  regions.  Upper  middle:  after  the  etching  of  the 
silicon  in  the  gate  region.  Lower  middle:  after  the  polysilicon  has 
been  deposited  and  defined.  Bottom:  The  completed  test  transistor. 
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Figu  re  1  -22:  Microphotograph  of  experimental  test  structures. 


Figure  1  -23:  SEM  cross-section  of  isolation  structure  with  the  largest 
fabricated  interface  radius,  r-  1.00  pm. 


the  substrate  was  implanted  with  boron  to  increase  the  surface  doping  to  1.3  x  10 16  cm'3.  They 
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Figure  1-24:  SEM  cross-section  of  isolation  structure  with  the  smallest 
fabricated  interface  radius,  r  =  0.57  n m. 


were  measured  with  VQS  =  0.05  V  and  Vsug  =  0.0  V.  Curve  A  is  from  a  transistor  whose  source/drain 
junctions  have  been  driven  in  below  the  oxide  corner  and  thus  obeys  simple  one  dimensional  MOS 
equations.  Its  measured  subthreshold  slope  is  0.56  V/dec.  Curve  B  corresponds  to  the  structure  of 
Figure  1-23  which  has  the  largest  corner  radius  of  the  fabricated  devices.  It  displays  a  slightly  flatter 
subthreshold  slope,  0.74  V/dec,  and  higher  threshold  voltage  compared  to  A.  The  transistor  whose 
trench  had  been  etched  with  Freon  115  and  had  all  the  field  oxide  thermally  grown,  has  the 
subthreshold  characteristics  of  curve  C.  Its  interface  radius  is  smaller  and  the  resulting  subthreshold 
curve  is  flatter  than  the  previous  two  curves,  0.80  V/dec.  The  fourth  subthreshold  curve,  labeled  D, 
corresponds  to  the  structure  of  Figure  1  -24.  This  device  has  the  sharpest  corner  of  all  the  fabricated 
devices.  As  a  result  the  subthreshold  curve  is  much  flatter  than  the  previous  curves  and  the  isolation 
is  significantly  improved.  The  subtheshold  slope  of  Curve  D  is  1.10  V/dec.  Between  devices  A  and  D, 
there  is  about  a  factor  of  two  improvement  in  isolation  where  isolation  is  defined  as  the  gate  voltage 
needed  to  produce  1  pA/pm  of  channel  width  of  leakage  current. 

The  effects  of  a  substrate  bias  are  shown  in  Fig  re  1  -26.  The  log(/Q)  vs.  VQS  curves  of  two  different 
devices  are  shown,  one  with  the  corner  effect  and  one  without.  The  structures  were  fabricated  on  a 
boron  implanted  substrate  and  have  channel  lengths  of  2.5  /*m  and  a  widths  of  500  fim,  VQS  =  0.05 
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Figu  re  1-25:  Plot  of  measured  I--V  characteristics  of  four  devices  with 

varying  amounts  of  the  comer  effect.  Curve  A  is  from  a  device  with 
no  corner  effect.  Curves  B,  C,  and  D  are  from  devices  with  decreasing 


interface  radii. 


.  =  0.0  V,  Vn-  =  0.05  V,  channel  length 


is  7.5  fim,  width  is  500  jim.  Nsuge&  1.3  X 


V.  Curves  A  and  B  are  from  a  transistor  without  the  corner  effect,  and  a  substrate  bias  of  0.0  and  -0.5 
V  respectively.  Curves  C  and  0  are  from  a  device  with  the  smallest  fabricated  corner  radius  (Figure 
1-24),  and  a  substrate  bias  of  0.0  and  -0.5  V  respectively.  The  shift  between  curves  A  and  B  is  3.4  V  at 
a  current  level  of  500  pA,  while  the  shift  between  curves  C  and  0  is  6.0  V.  The  characteristics  of  the 
comer  region  produce  a  larger  shift  in  the  threshold  voltage  for  a  given  substrate  bias,  and  hence 
improved  isolation.  The  subthreshold  slopes  also  steepen  for  both  cases  with  the  addition  of  the  bias, 
as  expected. 


The  SEMs  of  Figures  1-23  and  1-24  can  be  used  to  create  GEMINI  simulations  with  the  same 
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Figure  1  -26:  Plot  of  measured  l-V  characteristics  of  two  devices 

displaying  the  effect  of  substrate  bias.  Curve  A:  no  comer  effect, 
vsub  1 0  0  Curve  B:  V$UB  »  -0.5  V.  Curve  C:  with  comer  effect, 

VS(jg  =  0.0  V;  Curve  D:  V$UB  ®  -0.5  V.  Channel  length  is 
2.5  pm,  width  is  500  pm.  There  is  a  larger  shift  in 
the  curves  due  to  applied  substrate  bias  in  the  device  with  the  corner 
effect.  NS(jgp*i  1.3  x  TO^cm"3. 

physical  structure.  The  structure  can  be  simulated  over  a  range  of  gate  voltages  to  create  a 
subthreshold  curve.  The  simulated  curves  agree  with  the  experimentally  measured  curves,  and  as 
can  be  seen  in  Figure  1-20,  the  GEMINI  simulations  agree  with  the  analytic  model.  It  should  be  noted 
that  it  is  difficult  to  ascertain  the  exact  value  of  the  interface  radius  from  the  SEMs  and  GEMINI 
simulations,  which  hinders  efforts  to  make  exact  comparisons. 

Figure  1-27  displays  the  insensitivity  of  a  transistor  with  the  corner  effect  to  short  channel  effects. 
The  subthreshold  curves  of  two  devices,  with  and  without  the  comer  effect,  are  shown  for  drain 
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voltages  of  0.5  and  5.0  V.  The  structures  have  channel  lengths  of  2.5  pm  and  widths  of  500  pm.  The 
substrate  is  uniformly  doped  with  a  resistivity  of  6  Q-cm,  and  the  substrate  bias  is  0.0  V.  Because  of 
the  more  lightly  doped  substrate,  the  threshold  voltages  are  much  smaller  than  in  the  previous  two 
figures.  Curves  A  and  B  are  for  a  device  with  deep  junctions  that  obeys  the  ID  MOS  equations  with 
drain  voltages  of  0.5  and  5.0  V  respectively.  Curves  C  and  D  are  for  a  device  with  the  corner  effect 
(Figure  1*24),  and  drain  voltages  of  0.5  and  5.0  V  respectively.  There  is  a  much  larger  shift  with  the 
increased  drain  voltage  for  the  normal  device  than  for  the  device  with  the  corner  effect.  This  is  due  to 
the  isolation  of  the  comer  effect  transistor's  potential  barrier  from  the  drain  voltage.  The  shifting  of 
the  barrier  from  the  center  of  the  channel  to  the  source  makes  drain  induced  barrier  lowering  (DIBL) 
much  less  important,  as  was  explained  earlier  and  illustrated  in  Figure  1-15. 

To  further  illustrate  die  comer  effect's  immunity  to  short  channel  effects,  consider  Figure  1-28.  The 
drain  current  is  plotted  versus  the  drain  voltage  for  the  same  two  devices  that  are  used  in  Figure  1-27. 
The  gate  voltages  are  held  constant  at  values  that  put  each  transistor  in  the  subthreshold  regime.  For 
an  ideal  long  channel  ID  device,  theory  predicts  /Q  proportional  to  1-exp {-qV0S/kT).  The  curve  for 
the  structure  without  the  comer  effect,  A;  shows  this  rise  in  current  for  small  values  of  VD.  However, 
as  the  drain  voltage  continues  to  increase,  the  current  begins  to  rise  at  an  increasing  rate.  This  is 
indicative  of  DIBL  as  the  increasing  VQ  begins  to  reduce  the  potential  barrier  height.  Curve  B,  of  the 
structure  with  the  comer  effect,  does  not  show  this  superlinear  increase  in  current,  but  instead  l0 
levels  off  as  is  increased.  Clearly,  the  comer  effect  can  significantly  improve  the  isolation  provided 
by  parasistic  transistors  with  small  channel  lengths. 

1 .2.1 .3  Summary  of  Results 

The  availability  today  of  controlled  anisotropic  etching  makes  it  feasible  to  consider  trench-like 
structures  where  the  isolation  oxide  extends  significantly  deeper  into  the  silicon  than  the 
source/drain  junctions.  When  this  is  done,  a  substantial  improvement  in  isolation  is  obtained  due  to 
geometrical  arguments.  The  corner  effect  provides  flatter  subthreshold  slopes  and  higher  threshold 
voltages,  and  is  strongly  dependent  on  the  junction  depth  and  the  interface  radius  of  curvature.  A 
mathematical  analysis  of  the  physics  of  the  corner  effect  results  in  an  equation  relating  the 
subthreshold  slope  and  the  radius  of  the  interface.  That  strong  dependence  of  the  corner  effect  on 
the  radius  was  verified  by  the  testing  of  several  test  transistors  with  different  radii.  These 
measurements  also  demonstrated  the  improved  isolation  properties  with  substrate  bias  and  drain 
voltage. 

With  this  better  understanding  of  the  device  physics  of  isolation  transistors,  some  observations  can 
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Figure  1  -27:  Plot  of  measured  l-  V  characteristics  of  two  devices 

displaying  the  effect  of  drain  voltage.  Curve  A:  no  corner  effect, 
vos  *  0  5  V;  Curve  S:  VQS  =  5.0  V.  Curve  C:  with  comer  effect, 

VDS  =  0.5  V;  Curve  0:  V ^  *  5.0  V.  Channel  length  is  2.5  /tm, 
width  is  500  jim.  There  is  a  smaller  shift  in 
the  curves  due  to  the  applied  drain  voltage  in  the  device  with  the 
comer  effect.  NS(JB=  1.5  x  J0,5cm'3. 

be  made  about  the  various  new  fabrication  techniques.  Isolation  that  is  made  by  etching  away  the 
Field  oxide  in  the  active  device  regions  [14],  or  by  selectively  oxidizing  a  deposited  poiysilicon  layer 
into  the  field  oxide  [10],  will  not  demonstrate  the  two-dimensional  comer  effect  described  in  this 
paper  because  the  Si/Si02  interface  is  planar  and  the  oxide  is  above  the  source/drains.  In  addition 
to  not  realizing  the  better  isolation  due  to  the  corner,  short  channel  effects  will  remain  a  major 
problem.  Other  isolation  techniques  that  bury  the  oxide  in  a  trench-like  structure,  can  take  advantage 
of  the  better  isolation  available  to  it  by  adjusting  the  relative  positions  of  the  oxide  interface  and  the 
source/drain  junction  depths.  In  addition  to  improved  immunity  to  short  channel  effects,  these 
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Figu  re  1  *28:  Plot  of  the  drain  current  versus  the  drain  voltage  with  a 

constant  gate  voltage.  Curve  A  is  from  a  device  with  no  comer  effect 
and  VQ  » -0.06  V.  Curve  B  is  from  a  device  with  the  comer  effect  and 
VQ  *  1.07  V.  Channel  length  is  2.5  pm,  width  is  500  pm.  The 
device  with  the  comer  effect  does  not  experience 
drain  induced  barrier  lowering.  Nc.,a  »  7.5  X 

to'W.  SU8 

techniques  provide  a  planar  top  surface  for  easier  subsequent  processing.  However,  in  addition  to 
optimizing  electrical  characteristics  through  a  careful  choice  of  isolation  region  geometry,  other 
effects  such  as  stress,  defect  production  in  the  substrate,  and  charges  on  the  sidewalls  should  also 
be  considered  [3].  This  paper  has  attempted  to  point  out  that  with  these  new  techniques,  both  device 
physics  and  processing  techniques  can  point  toward  the  same  goal.  Burying  the  oxide  and 
sharpened  oxide  comera  means  improved  isolation  and  more  compact  VLSI  layouts. 
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1 .2.2  Alternative  Gate  Dielectrics  for  Submicron  MOS  VLSI 
Investigators:  M.  M.  Moslehi  and  K.  C.  Saraswat 

1.2. 2.1  Scientific  Objective 

Due  to  the  continuing  growth  in  the  integration  density  of  integrated  circuits  and  the  technological 
and  reliability  problems  with  silicon  dioxide  in  the  very  thin  regime  (  <  100  angstroms),  there  is  a 
demand  for  new  higher  quality  ultrathin  insulators  to  replace  silicon  dioxide  in  highly  scaled  down 
VLSI  IGFETs  and  memory  devices  such  as  DRAM  and  EEPROM.  Thermal  nitridation  of  Si  and  thin 
SiOs  has  been  investigated  in  the  past  few  years  and  appears  to  be  an  alternative  to  the  oxidation 
process  to  grow  good  quality  films  in  the  very  thin  regime  ( <  100  angstrom). 

Most  of  the  recent  work  on  thermal  nitridation  has  been  on  technological  issues  such  as 
composition,  growth  kinetics,  oxidation  resistance,  and  interactions  with  impurity 
diffusion  [1, 17, 14, 4, 3, 11]. 

In  the  last  report  of  this  program,  we  reported  a  new  technique  to  thermally  nitride  Si  and  S/02  in  a 
cold-wall  rf-heated  reactor  in  an  ammonia  gas  ambient.  Growth  kinetics  of  the  nitride  and  nitroxide 
films  was  studied  in  the  temperature  range  of  950  to  1200°C  [14].  The  physical  properties  of  these 
films  were  characterized  by  auger  electron  spectrometry,  Rutherford  back  scattering,  and  etch  rate 
measurements  [14].  Some  electrical  characterization  was  done  by  the  capacitance- voltage 
measurement  technique  [13].  In  a  few  cases  devices  made  with  nitroxide  and  nitride  have  been 
electrically  characterized  [7, 8, 13].  However,  most  of  the  work  in  the  past  has  been  on  nitrided  oxide. 
Nitridation  of  oxide  results  in  die  exchange  reaction  of  nitridation  species  with  oxide  at  the  surface 
and  interface  regions  and  to  a  lesser  extent  in  the  bulk  region.  It  has  been  reported  that  the 
generation  of  radiation-induced  interface  states  in  MOS  devices  is  eliminated  after  proper  high 
temperature  nitridation  of  oxide  [15].  This  may  be  attributed  to  structural  changes  in  the  strained 
region  of  the  oxide  after  nitridation  in  ammonia.  We  have  observed  [13],  as  ito  et  a!  [9]  have  also 
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reported,  that  the  breakdown  characteristics  of  MOS  structures  can  be  improved  by  annealing  the 
S/02  films  in  ammonia  gas  at  high  temperatures.  The  effects  of  nitridation  of  S/02  in  ammonia  on 
electron  trapping  has  been  investigated  [12, 2]  and  it  has  been  reported  that  ammonia  nitridation 
results  in  a  significant  increase  in  electron  trapping.  However,  the  data  presented  so  far  have  not 
been  sufficient  to  draw  firm  conclusions  regarding  trapping  in  nitroxide  insulators.  Moreover,  effects 
of  annealing  of  nitroxide  insulators  in  nitrogen,  oxygen,  or  argon  ambients  on  trapping  characteristics 
should  be  studied  more  in  detail. 

Thermal  nitride  films  can  be  grown  by  high  temperature  nitridation  of  silicon  in  pure  ammonia  [14], 
ammonia  plasma  [16, 10],  or  nitrogen-hydrogen  plasma  [6].  Ito  et.  ai.  have  fabricated  IGFETs  with 
thermal  nitride  gate  insulators  showing  good  breakdown  characteristics  and  high 
transconductance  [7, 10, 6].  However,  a  more  detailed  electrical  characterization  of  thermal  nitride 
films  including  trapping  has  not  been  reported  in  the  past.  In  this  report  the  results  of  a  study  on 
conduction  and  trapping  in  thermal  nitride  gate  insulators  is  presented.  Cross-sectional  TEM  has 
been  used  to  study  the  structure  of  thermal  nitrides  and  their  interface  with  the  silicon  substrate. 

1. 2.2.2  Progress 

Experimental  Procedure 

MIS  devices  with  thin  oxide,  nitroxide,  and  nitride  gate  insulators  were  fabricated  with  aluminum 
gate  electrodes  on  n-  and  p-type  Si  substrates.  The  S/02  used  in  this  work  was  grown  to  a  thickness 
of  about  100  angstroms  in  a  furnace  in  oxygen  at  900  °C.  All  thermal  nitridations  were  performed  in 
pure  ammonia  in  a  cold- wall  RF-heated  reactor  described  elsewhere  [13],  Nitridations  of  Si  and  S/02 
were  performed  at  temperatures  from  950  °C  up  to  1230  °C  for  times  ranging  from  30  minutes  to  4 
hours.  The  gate  areas  were  2.01  x  10'2,  5.03  X  10"3, 1.26  x  10*3,  3.14  x  10-4,  and  7.85  X  10*5  cm2. 
The  gate  electrode  was  Al-(1  %)Si  deposited  in  a  flash  evaporator. 

Results  and  Discussion 

Figure  1-29  illustrates  the  conduction  characteristics  of  an  MIS  device  with  64  angstroms  of  thermal 
nitride  grown  at  1230  °C  for  4  hours  on  n-type  (100)  silicon  with  a  resistivity  of  0.5  ohm.cm.  at  1230  °C 
for  4  hours  and  a  gate  area  of  7.85  x  10  5  cm2.  Fowler- Nordheim  (F-N)  behavior  can  be  observed 
over  about  six  decades  of  current.  At  very  high  current  levels  bulk  effects  on  conductivity  become 
important  resulting  in  deviation  from  Fowler- Nordheim  behavior.  The  F-N  parameters  can  be  obtained 
from  a  comparison  of  the  experimental  data  with  the  F-N  expression  J  =  a.E 2  exp(  b  /  E). 

It  is  well  known  that  for  thick  Si02  gate  insulators  the  current  density  at  a  given  field  is  independent 
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1/V(volt) 

Figu  re  1  -29:  l-V  characteristic  of  an  MIS  device  with  64  angstrom  nitride 

of  oxide  thickness.  However,  this  is  not  true  for  thermal  nitrides  in  the  ultrathin  regime  when  direct 
tunneling  through  the  nitride  becomes  important  Figure  1  -30  shows  the  l-V  characteristics  of  Al-gate 
devices  (e-beam  deposited)  with  gate  areas  of  7.85  X10'5,  3.14  x  10'4,  and  1.26  x  10  3  cm2.  The 
gate  nitride  was  grown  at  1100  °C  for  2  hours  resulting  in  about  40  angstroms  of  thermal  nitride.  By 
comparing  this  figure  with  that  for  64  angstrom  nitrides,  we  observe  a  sharp  increase  in  conductivity 
through  the  nitride  by  going  from  64  angstroms  to  40  angstroms.  The  increase  in  conductivity  is  too 
high  to  be  just  attributed  to  the  effect  of  a  reduction  in  thickness  of  the  nitride  on  F-N  tunneling.  In  fact 
for  nitride  films  of  40  angstroms  or  thinner  the  effect  of  direct  tunneling  on  conductivity  becomes  a 
dominant  factor.  The  excess  amorphous  silicon  content  of  thinner  nitride  films  might  be  an  additional 
factor  helping  to  increase  the  conductivity.  The  current  densities  in  devices  with  three  different  areas 
are  the  same  which  is  an  indication  of  no  edge  effects  on  conductivity  even  in  devices  with  gate 
nitride  as  thin  as  40  angstroms.  The  effect  of  direct  tunneling  in  thermal  nitride  films  thinner  than  40 
angstroms  implies  that  nitride  films  thinner  than  40  angstroms  may  result  in  high  gate  leakage 
currents  in  IGFETs  and  low  retention  in  EEPROMs.  This  was  also  verified  by  testing  the  devices  with 
nitride  insulator  thicknesses  ranging  from  40  angstroms  down  to  20  angstroms. 

Figure  1-31  shows  the  conduction  characteristics  of  devices  with  nitroxide  and  Si02  gate  insulators 
of  gate  electrode  area  5.03  X  10'3  cm2. 100  angstroms  Si02  was  grown  in  dry  02  at  900  °C  and  then 
nitrided  at  1000  °C  for  4  hours.  It  can  be  observed  that  as  a  result  of  thermal  nitridation  there  is  an 
increase  in  conductivity  through  the  device.  Assuming  that  the  current  is  mainly  due  to  injection  of 
electrons  from  the  gate  electrode  (device  in  accumulation  region),  the  nitridation  enhanced 
conductivity  can  be  attributed  to  lowering  of  the  barrier  height  at  the  surface  due  to  the  exchange  of 
oxygen  with  nitrogen  resulting  in  an  enhanced  tunneling  probability. 
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Figu  re  1  -30:  I- V  characteristics  of  devices  with  40  angstrom  nitride 

-6 


V(vo|ts) 

Figure  1*31:  Conduction  characteristics  of  devices  with  100  angstrom 

oxide  and  nitroxide 


Figure  1 -32(a)  illustrates  a  high  frequency  C-V  plot  of  an  MOS  device  with  100  angstroms  S/0? 
before  and  after  current  injection  of  1  /iA  for  15  seconds  through  the  insulator  (  device  in 
accumulation).  The  C-V  curve  was  shifted  to  the  right  by  about  0.5  volt  which  is  an  indication  of 
electron  trapping.  Figure  1  -32(b)  shows  the  high  frequency  C-V  characteristic  of  a  device  with  about 
60  angstroms  thermal  nitride.  The  C-V  curve  exhibited  no  hysteresis  or  shift  during  measurement  and 
also  after  comparable  charge  injection  indicating  extremely  low  carrier  trapping. 
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even  nitroxide.  Figure  1-33  shows  the  breakdown  characteristic  of  devices  with  nitride  gate  insulators 
on  n-type  silicon.  No  destructive  breakdown  could  be  observed  for  currents  up  to  a  few  hundred  pA. 
The  breakdown  measurements  on  a  device  were  reproducible  indicating  the  nondestructive  nature  of 
breakdown  for  even  very  high  current  densities.  This  will  be  more  clear  when  the  time  dependent 
electrical  breakdown  results  are  explained.  Under  the  same  conditions  devices  with  100  angstroms 
oxide  and  nitroxide  showed  destructive  breakdown  for  much  lower  current  densities.  Nitroxide  films 
were  observed  to  withstand  higher  current  densities  than  oxides  before  destructive  breakdown.  The 
breakdown  field  for  thermal  nitrides  varies  from  about  10  MV/cm  up  to  about  20  MV/cm  depending 
on  the  electrode  material,  device  area,  and  nitridation  conditions. 


Figu  re  1  -33:  Breakdown  characteristics  of  devices  with  thermal  nitride  insulators 

Conduction  characteristics  can  be  influenced  by  trapping.  The  F-N  plots  in  Figure  1  -34  illustrate  the 
effect  of  trapping  on  conduction  in  devices  with  60  angstroms  thermal  nitride.  I-V  measurements 
were  taken  before  and  after  injecting  25  juA  into  the  gate  for  100,  200,  and  500  seconds.  The 
rightmost  curve  shows  the  original  conduction  characteristic.  Upon  charge  injection  the  l-V  curves 
were  originally  shifted  slightly  to  the  left  and  then  again  back  to  the  right.  The  shifts  are  very  small  and 
no  destructive  breakdown  was  observed  even  after  500  seconds  current  injection  corresponding  to 
about  160  C/cm2  at  a  current  density  of  0.32  A/cm2.  This  is  much  larger  than  the  breakdown  charge 
density  in  100  angstrom  thermal  oxides  at  comparable  current  densities  which  was  observed  to  be 
about  20  C/cm2. 

Time  dependent  electrical  breakdown  characteristics  of  devices  with  thermal  nitride  gate  insulators 
were  also  examined.  Figures  1 -35(a)  and  1 -35(b)  show  these  characteristics  for  devices  with  nitride 
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gate  gate  insulators  grown  at  1230  °C  for  1  hour  and  4  hours  on  p-type  and  n-type  silicon.  For  an 
injected  current  density  of  1.27  C/cm2,  the  devices  with  nitride  grown  for  1  hour  on  p-type  silicon 
showed  destructive  breakdown  after  400  seconds  corresponding  to  about  510  C/cm2  whereas  the 
devices  with  nitride  grown  for  4  hours  did  not  show  any  destructive  breakdown  even  after  1000 
seconds  or  injection  of  1270  C/cm2.  These  charge  densities  are  much  higher  than  that  for  100 
angstrom  thermal  oxides  at  comparable  stress  current  densities  which  is  less  than  20  C/cm2.  In  both 
figures  the  insets  show  the  expanded  voltage  scale  plots  indicating  hole  trapping  in  the  first  case  and 
electron  trapping  in  the  second  one.  The  rate  of  trapping  is  extremely  low  which  can  be  attributed  to 
the  high  quality  of  the  thermal  nitride  gate  insulators  and  importance  of  field  emission  tunneling 
detrapping  mechanisms  in  ultrathin  gate  insulators. 

The  interface  quality  and  roughness  has  a  direct  relationship  with  the  surface  mobility  of  carriers  in 
MISFETs.  The  interface  roughness  observed  for  thermal  nitride  is  as  good  as  those  observed 
previously  [5]  for  Si/SiO^  structures.  We  measured  a  thickness  of  41  angstroms  for  this  film  using  an 
ellipsometer  with  the  index  of  refraction  fixed  at  2.0.  Using  the  cross-sectional  TEM  image  of  the  same 
film  and  the  (111)  interplanar  spacing  as  the  measurement  unit,  the  thickness  is  measured  to  be  43  =F 
1.6  angstroms  which  is  in  good  agreement  with  the  ellipsometry  result. 

1.2. 2. 3  Summary 

In  summary,  we  have  studied  the  electrical  characteristics  and  interface  morphology  of  thermal 
silicon  nitride  films.  The  results  indicate  extremely  low  trapping  in  these  films.  Films  thicker  than  40 
angstroms  can  be  used  as  gate  insulators  of  field  effect  transistors;  however,  those  thinner  than  that 
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Fl9ur.1-35:  Tim,  dependent  electrical  breakdown  characteristics  lor  devices  with 
nitnde  grown  at  1230  °C  tor  (a)  1  hour  on  p-type  and  (b)  4  hours 

on  n-type  Si 

show  high  leakage  doe  to  the  quantum  mechanical  direct  tunneling  limitation.  The  Interface 
abruptness  and  quality  is  comparable  to  high  quality  dry  Si02.  High  density  EEPROM  and  DRAM 

circurts  as  wet,  as  scaled  VLS,  IGFETs  should  find  thermal  silicon  nitrides  very  useful  tunneling  and 
gate  insulators. 
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1.2.3  Electronic  Properties  of  Polycrystalline  Silicon 

Investigators:  O.B.  Kao  and  K.C.  Saraswat 

1. 2.3.1  Scientific  Objective 

The  major  objective  of  this  study  is  to  investigate  the  oxide  fixed  charge,  Qp  and  the  interface 
traps,  Ng,  in  polysilicon  gate  MOS  devices  and  specifically  to  study  the  density  of  QF  and  as  a 
function  of  high  temperature  treatments  in  various  ambients  and  the  effect  of  device  dimension 
scaling. 

This  work  is  important  because  most  MOS  devices  today  have  Poly-Si  gates,  while  most  of  the  work 
on  Qf  annealing  has  been  based  on  data  from  aluminum  gate  devices  where  the  bare  gate  oxide  is 
quite  sensitive  to  processing  conditions  before  the  final  low  temperature  metalization  step.  The  QF 
triangle,  as  shown  by  Deal.et  a!.  [2]  illustrates  that  these  charges  vary  with  processing  temperature 
and  ambient.  In  an  earlier  experiment  [3],  we  reported  that  large  area  Poly-Si  gate  capacitors  are 
insensitive  to  processing  conditions,  because  polysiiicon  shields  the  gate  oxide  from  02  in  the 
ambient.  That  finding  led  us  to  pursue  the  edge  effects  of  scaled  down  Poiy-Si  gates,  especially  the 
peripheral  0F  generated  by  the  diffusion  of  oxygen  at  the  edges  of  the  gate.  This  peripheral  QF  will  be 
important  for  total  Qf  control  in  small  devices. 

Our  results  indicate  that,  even  for  gate  width  as  small  as  1  25  pm,  QF  underneath  the  polysilicon 
gate  is  unaffected  by  further  processing  steps,  including  high  temperature  oxidizing  ambients.  A 
two-dimensional  oxygen  diffusion  model  is  proposed  to  explain  this  phenomenon. 
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1. 2.3.2  Progress 

Experimental 

Polysilicon  gate  MOS  capacitors  with  finger  structures  were  fabricated  using  a  duai  oxide  approach 
with  a  gate  oxide  of  630  angstroms,  and  thick  field  oxide  of  7400  angstroms.  Ail  capacitors  have  the 
same  area.  The  width  of  the  fingers  and  the  spacing  between  the  fingers  were  kept  equal  and  varied 
from  10  jim  to  1.25  pm  resulting  in  a  larger  perimeter  for  smaller  finger  widths.  The  wafers  were 
annealed  in  02  and  Ar  at  different  high  temperatures  and  for  different  lengths  of  time.  After  each  high 
temperature  anneal,  the  wafers  were  annealed  in  forming  gas  at  450°C  for  30  minutes,  before  C-V 
measurements  were  taken  to  determine  VFB  and  Qp.  From  the  difference  between  high-frequency 
and  quasi-static  data,  Dn  was  determined  to  be  less  than  2x1 010  cm’2.  The  low  Dn  value  indicates  that 
hydrogen  in  the  forming  gas  diffuses  through  polysilicon  grain  boundaries  and  anneals  the  interface 
traps  sufficiently.  Figure  1-36  shows  a  series  of  high  frequency  C-V  curves  for  devices  with  different 
perimeter  to  area  ratios.  The  fringing  capacitance  is  evident  when  the  devices  are  in  inversion, 
resulting  in  the  different  values  for  Cm/r).  It  is  negligible  in  accumulation.  This  is  understandable 
because  in  accumulation,  the  oxide  capacitance  dominates  and  the  fringing  capacitance  is  negligible 
because  of  the  thin  gate  oxide.  When  the  device  is  in  inversion,  however,  silicon  depletion 
capacitance  dominates.  In  this  case,  the  fringing  field  extends  the  depletion  region  outside  the  gate 
and  reduces  the  maximum  depeletion  width,  thus  increasing  the  silicon  depletion  capacitance  and 
the  total  capacitance.  The  inversion  region  fringing  capacitance  is  proportional  to  the  perimeter  to 
area  ratio. 

Results  and  Discussion 

Experimental  results  of  annealing  are  shown  in  Figure  1-37.  The  normalized  VFB(VFfl/VFBftn(Ha())  is 
plotted  for  a  sequence  of  successive  annealing  conditions.  It  can  be  seen  that  the  capacitors  do  not 
show  any  significant  changes  in  VFfl  or  QF  after  anneal  in  02  at  900°C  and  1000°C  for  different 
periods  of  time  and  different  02  partial  pressures,  or  after  a  final  Ar  anneal  at  1000°C.  This  is  true  for 
devices  of  all  dimensions,  including  the  ones  with  1.25  pm  wide  fingers.  This  result  indicates  that 
there  is  negligible  change  of  0F  in  the  gate  oxide  underneath  the  polysilicon  gate.  Since  0F  is  a 
product  of  the  oxidation  process,  a  plausible  explanation  is  that  there  is  simply  very  little  oxygen 
underneath  the  polysilicon  gate,  which  may  be  as  narrow  as  1.25  pm.  This  seems  to  be  in  obvious 
contradiction  to  what  would  be  predicated  by  the  high  diffusivity  of  02  in  S/02,  about  1CT8  cm2/sec. 
To  resolve  this  dilemma,  a  study  of  the  diffusion  of  02  in  the  gate  oxide  underneath  polysilicon  was 
undertaken. 
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Figure  1  -36:  High-frequence  C-V  plots  for  different  finger  widths  (L), 
or  perimeter  to  area  ratiois  (P/A).  The  inversion  region 
capacitance  is  proportional  to  P/A. 
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Figu re  1  -37:  Normalized  flat  band  voltages,  '<VFB(norm)>,  defined  as 

KVpg/V FB(initiai)^ for  a  series  of  annealing  conditions.  There  is 
negligible  change  in  caused  by  peripheral  QF  under  Poly-Si 

gates. 


The  Deal-Grove  oxidation  model  [2]  can  be  extended  to  two  dimensions  [1].  Our  case  is  a  two- 
dimensional  steady-state  diffusion  problem  with  a  source  and  different  sinks: 
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DV*C  =  =  0  Laplace' »  Eq. 

ot 

D  —  =  h(C *  -  (7)  at  Oj  tource 
dC 

D—=ktC  at  Si/SiOj  interface 
dC 

D—  =  0  at  non  —  oxidizing  interface 


Numerical  solution  was  carried  out  with  an  iterative  finite-difference  boundary  value  method.  Figure 
1-38  shows  the  simulated  contour  plots  of  02  concentration  near  the  edge  of  a  nitride  cap  (non¬ 
oxidizing  boundary)  and  of  a  polysilicon  gate  (oxidizing  boundary),  respectively.  Figure  1-39  shows 
the  02  concentrations  at  the  S/02/S/-substrate  interface,  for  both  the  nitride  case  and  polysilicon 
case.  It  can  be  seen  that  even  underneath  nitride,  the  steady  state  02  concentration  under  the  cap  is 
relatively  small,  despite  the  high  diffusivity  of  02,  consistent  with  result  shown  in  [1].  This  makes 
sense  because  the  oxidizing  boundary  of  the  Si  surface  serves  as  a  sink  for  02,  limiting  its  distribution 
deep  underneath  the  gate.  The  fast- oxidizing  polysilicon  consumes  even  more  02,  thus  pulling  most 
02  away  from  the  S/02/S/- substrate  interface,  leaving  even  less  02  available  for  Qf  generation  in  the 
interface. 

1. 2.3.3  Summary 

In  conclusion,  our  study  shows  that  the  polysilicon  gate  not  only  acts  as  an  oxygen  barrier  above 
the  gate  oxide,  it  also  keeps  oxygen  away  from  the  S/02/S/-substrate  interface  under  the  gate  edges. 
Even  for  very  small  polysilicon  gate  dimensions,  QF  anneals  under  the  gate  independent  of  previous 
process  history.  In  addition,  subsequent  high  temperature  exposures  to  oxidizing  ambients  have  no 
effect  on  Qp.  However,  as  devices  are  scaled  down  further,  other  edge  effects  may  become 
important:  such  as  fringing  fields,  and  non-uniform  gate  oxide  due  to  polysilicon  oxidation. 

1. 2.3.4  References 

1 .  Chin.D.  J.,  Oh,S.  Y.,  Hu,S.  M.,  Dutton, R.  W.,  and  Moll.J.  L.  "Two- Dimension  Oxidation  Modeling." 
IEEE  Trans,  on  Electron  Devices  ED-30  (July  1983). 

2.  Deal, B.E.  and  Grove, A.S.  "General  Relationship  for  the  Thermal  Oxidation  of  Silicon.”  J.  Appl. 
Phys.  (December  1965). 
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3.  McVittie.J.  P.,  Chang, G.,  Kao.O.  B.,  and  Patton.G.  L.  Annealing  of  Fixed  Charges  in  Poly-Si  Gate 
MOS.  ECS  Meeting  Extended  Abstract,  May,  1983. 

1. 2.3.5  Recent  JSEP  Publications 

1 .  Kao.D.B.,  Saraswat.K.C.,  and  McVittie.J.P.  Annealing  of  Oxide  Fixed  Charges  in  Scaled 
Polysilicon  Gate  MOS  Structures.  IEEE  Semiconductor  Interface  Specialist  Conference, 

December,  1983. 

2.  Kao.D.B.,  Saraswat,K.C.,  and  McVittie.J.P.  Annealing  of  Oxide  Fixed  Charges  in  Scaled 
Polysilicon  Gate  MOS  Structures,  submitted  to  IEEE  Transactions  on  Electron  Devices 

3.  McVittie.J.  P.,  Chang, G.,  Kao.D.  B.,  and  Patton.G.  L.  Annealing  of  Fixed  Charges  in  Poly-Si  Gate 
MOS.  ECS  Meeting  Extended  Abstract,  May,  1983. 

1.2.4  Ohmic  Semiconductor  Contacts  with  Low  Minority-Carrier  Recombination  Velocity 
Investigator:  R.  M.  Swanson 

1 .2.4.1  Scientific  Objectives 

This  research  program  has  aimed  at  developing  SIPOS  (Semi- Insulating  Poly  Silicon)  emitters  for 
bipolar  transistor  applications. 

.  -if-.  «.  1  ' 

4  y 

1. 2.4.2  Progress 

Using  technology  improvements  discussed  in  previous  reports  the  emitter  Gummel  number  was 
increased  two  orders  of  magnitude  over  that  attainable  with  a  conventional  diffused  emitter.  Due  to 
this  success,  alternative  support  has  now  been  obtained  to  begin  looking  at  the  impact  these  emitters 
might  have  the  performance  of  high  speed  bipolar  devices.  A  considerable  speed  increase  is 
expected  through  the  use  of  SIPOS  emitters. 

The  relaxation  of  design  constraints  that  result  from  the  reduction  of  J^  through  SIPOS  emitters 

provides  greater  flexibility  and  may  suggest  new  guidelines  for  device  optimization.  It  has  been 

demonstrated  that  shallow  homojunctions  made  in  silicon  are  limited  by  materials  parameters 

associated  with  high  level  doping  effects  to  values  in  excess  of  IxlO'12  A/cm2.  Examination  of 

parameters  on  a  phenomenlogical  level  indicate  that  several  alternate  structures  can  reduce  J„ 

00 

several  orders  of  magnitude  below  this  "fundamental"  limit;  however,  the  SIS  structure  incorporating 
a  thin  interfacial  silicon  dioxide  that  both  passivates  the  interface  and  imposes  a  transport  barrier  to 
holes  appears  to  offer  a  solution  that  is  compatable  with  standard  process  technology.  Modeling  of 
both  static  and  dynamic  behavior  of  this  structure  produced  guidelines  regarding  the  extent  of 
outdiffusion  that  could  be  tolerated.  It  was  seen  that  reduction  of  J^  involved  attaining  a  profile  with 


JSEP  Annual  Report 


Page  63 


June  1984 


Stanford  University 


Stanford  Electronics  Laboratories 


sufficiently  high  surface  concentration  to  achieve  a  low  equilibrium  surface  hole  density,  yet  shallow 
enough  to  minimize  the  impact  of  Auger  recombination.  High  frequency  performance  which  relies  on 
the  minimization  of  hole  storage,  was  shown  to  favor  the  shallow  structures  as  well.  Reduction  of  the 
extent  of  the  outdiffused  layer  consistent  with  the  surface  concentration  requirement  was  shown  to 
favor  both  static  and  dynamic  performance. 

Experimental  BJTs  incorporating  Poly  and  SIPOS  contacts  were  used  to  develop  the  fabrication 
technology  for  achieving  low  static  values  of  It  was  shown  that  for  POLY  emitters,  that  the 
presence  of  the  interfacia)  oxide  is  mandatory  for  good  performance.  "Intimate"  POLY  contacts 
behaving  essentially  as  extensions  of  the  monocrystalline  regions,  do  not  seem  to  be  able  to  offer 
much  in  the  way  of  achieving  low  absolute  values  of  JM  although  relative  enhancements  are  readily 
attained.  Although  no  attempt  was  made  to  understand  hole  transport  mechanisms  in  heavily  doped 
POLY,  it  is  felt  that  they  are  understood  albeit  on  an  unsophisticated  phenomenlogical  level.  The 
situation  for  SIPOS  contacts  is  different  Low  oxygen  content  films  behaved  surprisingly  well. 
Although  all  indications  point  to  an  interfacial  barrier  action  for  this  type  of  contact  as  well,  the  exact 
nature  and  even  origin  of  this  interfacial  layer  is  poorly  understood.  In  part,  this  is  attributable  to  a 
lack  of  adequate  control  over  the  SIPOS  deposition  ambient. 

1. 2.4.3  Future  Plans 

Further  investigation  may  serve  to  clarify  some  of  these  issues.  Since  the  most  interesting  physics 
appears  to  lie  in  the  low  oxygen  content  SIPOS  films,  tighter  control  over  residuals  in  the  deposition 
ambient  is  necessary.  A  more  sophisticated  wafer  loading  system  is  needed  to  limit  the  growth  of 
unintentional  exides,  and  special  care  must  also  be  taken  to  avoid  contamination  by  residual 
organics.  If  the  interfacial  oxide  is  found  to  result  from  the  annealing  of  the  SIPOS  film  itself,  then  the 
kinetics  of  this  reaction  need  to  be  studied.  The  long  term  stability  of  these  interfacial  layers  has  not 
been  established;  work  on  possible  failure  mechanisms  is  needed.  Characterization  of  majority 
carrier  series  resistance  was  only  cursorily  done,  and  only  for  the  SIPOS  contacts.  This  is  a  major 
drawback,  as  it  makes  evaluation  of  the  relative  merits  of  SIPOS  and  POLY  contact  structures 
difficult.  This  is  certainly  an  area  that  needs  extensive  characterization  before  application  of  these 
structures  to  high  speed  or  high  current  density  applications  can  be  considered  practical. 

The  dynamic  performance  of  these  contact  structures  was  not  experimentally  investigated.  Actual 
device  geometries  optimized  to  take  advantage  of  the  low  J  need  to  be  fabricated  to  see  if  the 
potential  gains  that  were  suggested  can  be  realized  in  practice.  For  application  to  the  point  contact 
solar  cell,  the  performance  of  the  P-type  counterpart  must  be  established.  Comparison  to  the  N-type 
may  also  shed  some  light  on  the  physics  of  these  contacts. 


JSEP  Annual  Report 


Page  64 


June  1984 


Stanford  University 


Stanford  Electronics  Laboratories 


1. 2.4.4  Recent  JSEP  Publications 

1.  Kwark.Y.  and  Swanson, R.M.  "Optica)  Absorption  of  Thin  SIPOS  Films. "  J.  Electrochemical  Soc. 

129. 1  (January  1982),  197-201. 

2.  Yablanovich,  E.,  Swanson,  R.  M.,  and  Kwark,  Y.  H.  An  n-SIPOS:p- SIPOS  Homojunctionn  and  a 
SIPOS-Si-SIPOS  double  heterostructure.  Proceedings  17th  IEEE  Photovoltaics  Specialists 
Conference,  1984. 

1 .2.5  Interaction  of  Transition  Metals  with  Silicon  Surfaces 
investigator:  W.  E.  Spicer,  I.  Lindau 

1. 2.5.1  Scientific  Objective 

During  this  reporting  period  we  have  completed  the  first  preliminary  investigation  of  the  interaction 
of  transition  metals  with  Si  surfaces.  The  primary  objective  was  to  study  interaction  on  well- 
characterized,  i.e.  atomically  clean  and  structurally  ordered,  Si  surfaces.  To  this  end  the  interfaces 
were  prepared  by  cleavage  of  Si  in  situ  with  subsequent  metal  deposition  under  ultrahigh  vacuum 
conditions. 

1 .2.5.2  Progress 

Mo  was  deposited  onto  the  in  situ  cleaved  Si(111)  2x1  surface  with  coverages  ranges  from  0.1  to 
20  atomic  layers  [1].  The  induced  changes  in  the  electronic  structure  was  studied  by  photoemission 
using  synchrotron  radiation  in  the  soft  x-ray  region.  This  technique  was  also  used  to  determine  the 
abruptness  of  the  interface.  Interestingly,  the  room  temperature  results  showed  that  the  interface  is 
abrupt  on  an  atomic  scale,  contrary  to  results  we  have  obtained  earlier  for  Pd/Si  and  Pt/Si  where 
intermixing  takes  place  over  several  layers.  In  the  latter  two  cases  a  small  amount  of  oxygen  (i.e. 
about  one  monolayer)  will  drastically  reduce  the  interdiffusion.  In  the  case  of  Mo  the  effect  of  oxygen 
is  much  less  pronounced  [2].  The  next  step  in  this  work  will  be  to  examine  the  Schottky  barrier 
heights  in  these  different  situations.  As  far  as  the  chemical  reactions  are  concerned  we  observe 
siiicide  formation  already  at  room  temperature  at  the  low  coverages.  (The  Mo  bulk  silicide  forms  at 
525®C.  Even  for  submonolayer  coverages  the  Mo-Si  bonding  has  the  same  electronic  characteristics 
as  the  bulk  silicide,  Mo-Si2.  Thus  on  an  atomically  local  level  there  is  a  trend  to  have  the  bulk  Mo-Si 
configuration. 
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1.2. 5.3  References 
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"Exploiting  Photon  Energy  Dependence  in  Photoemission  from  Si(t  1 1  )-Mo  Interface.”  Physica  B  795 
(1983),  117-118. 

2.  Rossi,  G.,  Abbati,  I.,  Lindau,  I.,  Nogami,  J.,  Spicer,  W.  E.  The  Electronic  Structure  of  Clean  and 
Oxygen  Exposed  Si(111)-Mo  Interfaces.  In  preparation  for  submission  to  Phys.  Rev.  B. 

1. 2.5.4  Recent  JSEP  Publications 

1 .  Rossi,  G.  and  Lindau,  I.  "Compound  Formation  and  Bonding  Configuration  at  the  Si/Cu 
Interface.”  Phys.  Rev.  B  28  (1983),  3597. 

1.2.6  Interaction  of  Impurities  with  Silicon  Surfaces  and  Interfaces 
Investigator:  C.  R.  Helms 

1. 2.6.1  Scientific  Objective 

In  past  periods  of  this  contract,  the  objective  of  this  part  of  the  program  was  to  determine  the  effect 
of  arsenic,  boron,  and  other  impurities  on  the  electronic  structure  of  the  silicon  free  surface.  The  first 
phase  of  this  study  was  completed  last  year  with  the  Ph.D.  Dissertation  of  S.  C.  Perino  entitled  "The 
Electronic  Structure  and  Deposition  Kinetics  of  Arsenic  on  the  Silicon  Surface"  [2].  In  this  period  we 
have  continued  our  studies  specifically  to  determine  the  effect  of  interfacial  contamination  layers  on 
the  Schottky  barrier  height  of  silicon-metal  contacts. 

1. 2.6.2  Progress 

In  this  period,  we  have  invetigated  the  effect  of  the  presence  of  thin  (<  10  angstrom)  interfacial 
layers  of  oxide  and  carbon  on  the  Schottky  barrier  height  for  selected  silicon  metal  contacts.  These 
contamination  layers  can  be  deposited  and  indeed  represent  residues  from  some  of  the  standard  1C 
Chemical  Cleaning  procedures  [1].  Of  interest  here  was  the  effect  of  the  thin  oxide  left  on  the  surface 
after  the  "RCA"  clean  with  the  final  HF  etch  left  out  and  the  effect  of  thin  carbon  residues  left  after  the 
full  RCA  clean  including  the  HF  etch  final  step.  We  find  that  the  "RCA  oxide"  residue  has  a  large 
effect  on  Schottky  barrier  height  whereas  the  carbon  residue  has  very  little  effect.  For  example,  the 
thin  RCA  oxide  (10 )  can  increase  (decrease)  the  Schottky  barrier  of  Ti  on  p(n)  type  silicon  from  oco.6 
(0.5)  eV  to  0.9  (0.2)  eV.  Although  this  is  desirable  for  some  applications,  this  effect  is  not  thermally 
stable  for  temperatures  cc200  °C  and  can  lead  to  unreliable,  unreproducibie  barrier  values.  The  initial 
stages  of  this  investigation  are  complete  and  follow  on  fundinng  for  further  studies  based  o  these 
results  with  the  goal  of  establishing  a  complimentary  MESFET  technology  has  been  obtained  from  the 
Semiconductor  Research  Corporation. 
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1 .2.6.3  References 

1.  Meindl,  J.  0.  etal.  Annual  Progress  Report,  1  April  1982  -  31  March  1983.  Tech.  RepL  2,  Stanford 
Electronics  Laboratories,  Stanford  University,  1983. 

2.  Perino,S.C.  The  Electronic  Structure  and  Deposition  Kinetics  of  Arsenic  on  the  Silicon  Surface. 
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1.2.7  Interconnections  for  High  Density,  High  Performance  VLSI 
Investigators:  R.  F.  W.  Pease,  P.  de  la  Houssaye,  G.  L.  Yoffe 

1 .2.7.1  Scientific  Objectives 

One  result  of  the  increasing  scale  of  integration  is  that  the  packing  density  and  performance  of 
VLSI  circuits  is  at  present  limited  by  interconnects  and  not  by  intrinsic  speed  of  the  active  devices. 
The  aim  of  this  project  was  to  investigate  ways  in  which  the  limits  posed  by  interconnects  might  be 
overcome  so  that  advances  in  the  more  fundamental  properties  of  the  semiconductor  devices  will 
once  again  lead  to  advances  in  terms  of  improved  system  performance. 

Interconnects  can  limit  the  performance  of  VLSI  circuits  in  a  number  of  different  ways.  The  large 
fraction  of  space  occupied  by  interconnects  leads  to  long  transmission  paths  for  signals.  Thus  in  the 
most  basic  case,  in  which  die  interconnect  acts  as  a  lossless,  dispersionless  line,  there  will  be  a 
propagation  delay  set  by  the  parameters  of  the  line.  At  present  I.C.  metallization  paths  are  far  from 
acting  as  such  lines  and  the  delays  are  considerably  longer  than  those  of  a  line  propogating  signals  at 
a  speed  c/2.  With  the  advent  of  fine  line  lithography  the  number  of  squares  of  long  metal  runners  has 
increased  and  the  resistance  of  the  line  has  become  appreciable  in  that  long  lines  act  as  distributed 
RC  networks  (each  with  a  characterister  time  constant).  Shorter  lines  still  act  a  lumped  capacitor  that 
must  be  charged  or  discharged  in  active  devices  [1].  In  circuits  with  large  currents  being  drawn  the 
problems  of  electromigration  can  be  serious  and  voltage  drops  can  also  be  significant. 

Among  the  proposals  for  ameliorating  the  condition  imposed  by  todays  metallization  are: 

1.  Improve  the  metallurgy  to  improve  resistance  to  electromigration. . 

2.  Operate  at  low  temperature  to  improve  electrical  conductance. 

3.  Exploit  the  vertical  dimension  to  improve  the  conductance  per  square  and  also  reduce 
the  area  taken  up  by  interconnect  and  hence  reduce  interconnect  length. 
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1. 2.7.2  Progress 

The  first  approach  was  to  investigate  the  possibilites  of  improved  resistance  to  electromigration  of 
aluminum  alloy  runners  by  selective  annealing  with  a  scanning  electron  beam.  The  principle  in  to 
induce  all  grain  boundaries  to  be  perpendicular  to  the  direction  of  current  flow.  Previous  work  has 
shown  that  such  a  "bamboo"  structure  leads  to  greatly  improved  mean  time  to  failure  [3]  but  that  the 
structure  occured  naturally  only  under  restricted  conditions. 

The  results  obtained  showed  morphological  evidence  of  some  local  bamboo  structure  but  the 
practicality  of  achieving  such  a  structure  repeatedly  over  long  lengths  of  metal  ' over  a  variety  of 
underlying  films  as  occurs  in  an  integrated  circuit)  seemed  remote. 

However  the  ability  to  use  the  heating  effect  of  the  beam  to  fashion  the  metal  led  to  other 
interesting  possibilities  with  much  greater  potential  impact  on  the  interconnect  problem.  One  of  the 
main  difficulties  of  exploiting  the  vertical  dimension  to  achieve  many  layers  of  interconnect  is  that  of 
achieving  a  planar  structure  containing  arbitary  pattens  of  metal  and  insulator.  In  particular  the 
problems  of  filling  "posthole"  vias  with  metal  has  proved  particularly  awkard.  The  possibility  of  using 
the  directed  energy  of  a  laser  or  electron  beam  to  bring  about  controlled  local  melting  of  the  metal 
appears  quite  promising.  We  have  investigated  one  suggested  scheme  depicted  in  Figure  1-40. 

The  experiments  carried  out  so  far  have'  employed  lines  rather  than  individual  holes.  The  oxide 
thickness  is  0.5  pm  and  the  metal,  pure  aluminum,  was  deposited  by  evaporation  and  patterned 
lift-off.  Raster  scanning  with  a  20  kV,  3  jim  diameter,  electron  beams  of  20  /iA  current  resulted  in 
local  melting  and  planarization  was  accomplished  in  some  instances  (Figure  1-41).  It  may  be  that 
transient  large  area  heating,  such  as  can  be  obtained  with  a  pulsed  laser,  amy  be  more  suitable  and 
that  undue  heating  of  underlying  layers  can  still  be  avoided.  Thus  it  appears  that  transient  heat 
treatment  which  has  been  so  successful  in  fashioning  semiconductors  may  also  have  a  big  impact  in 
solving  some  of  the  metallization  problems.  Any  successful  realization  of  this  technology  will  require 
fundamental  studies  of  the  insulator  surface  and  of  the  metal/insulator  interface.  Nonetheless,  the 
technonlogical  rather  than  scientific  emphasis  of  this  project  as  it  is  now  configured  is  leadinng  us  to 
seek  support  from  programs  with  more  directed  aims  than  JSEP. 

1. 2.7.3  References 

1 .  Sinha,  A.  K.,  Cooper,  J.  A.,  Jr.,  and  Levinstein,  H.  J.  "Speed  Limitations  Due  to  Interconnect  Time 
Constants  in  VLSI  Integrated  Circuits."  IEEE  Elec.  Dev.  Lett.  EDL-3  (April  1982),  90-92. 
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deposited  metal  following  lift-off 


deposited  metal  following  local  melting 


Figu re  1  -40:  Steps  in  process  leading  to  a  planar  structure  containing 
arbitrary  patterns  of  metal  and  insulator.  The  critical  steps  are 
depositing  sufficient  metal  to  fill  etched  features  in  the  insulator 
and  controllably  melting  the  deposited  metal  to  planarize  the  top 
surface  without  damaging  underlying  layers  [2]. 
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Figure  1-41:  Scanning  electron  micrograph  of  aluminum  lines  in  1  pm-wide 
Si02  trenches  locally  melted  (in  lower  1  /3  of  picture)  to 
form  a  nearly  planar  structure 


2.  Tuckerman,  D.  B.  Unpublished  Work. 

3.  Vaidya,  S.,  Fraser,  D.  B.  and  Sinha,  A.  K.  Electromigration  Resistance  of  Fine-Line  Al  for  VLSI 
Applications.  Proceedings  18th  IEEE  Reliability  Physics  Symposium,  April,  1980,  pp.  165-170. 
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1.2.8  Study  of  Compound  Semiconductor  Device  and  Materials  Physics 

Investigator:  James  Harris 

1. 2.8.1  Scientific  Objective 

The  objective  of  this  project  is  to  investigate  heterojunction,  quantum  well  and  superlattice 
concepts  and  their  application  to  new  electronic  devices.  The  long  range  objective  is  to  study 
devices  utilizing  structures  of  ultra-smalt  dimensions  where  quantum  size  effects  totally  dominate  the 
transport  properties  and  energy  states  within  the  structures.  These  new  device  concepts  may  provide 
the  foundation  for  an  electronics  technology  with  vastly  superior  performance  compared  to  devices 
based  upon  current  device  concepts. 

1. 2.8.2  Progress 

Molecular  Beam  Epitaxy  (MBEJcan  readily  produce  structures  of  the  required  dimensions  (<200 
Angstroms)  and  is  thus  a  major  facility  required  to  carry  but  research  in  this  area.  Since  this  is  a  new 
project  at  Stanford,  the  first  year  objectives  were  to  define  a  program  in  lll-V  compound 
semiconductor  devices,  acquire  funding  for  an  MBE  system  and  begin  to  develop  the  necessary 
processing  and  characterization  facilities  to  carry  out  research  in  this  area. 

One  of  the  major  accomplishments  this  year  was  the  award  of  a  DARPA  program  to  fund  purchase 
of  an  MBE  system.  A  two  chamber  Varian  MBE  system  will  be  delivered  in  Aug.  1984.  This  system  will 
be  available  to  carry  out  the  research  on  ultra-small  dimension  structures  proposed  in  this  project. 

Examination  of  various  device  concepts  during  die  past  year  has  led  us  to  the  conclusion  that  the 
next  frontier  in  device  concepts  will  not  be  an  evolutionary  one,  but  a  revolutionary  one.  That  is,  even 
the  new  hetero-junction  concepts,  such  as  the  MODFET  and  heterojunction  bipolar  transistor,  are  still 
largely  based  upon  conventional  device  physics  and  carrier  transport  concepts.  These  structures  will 
certainly  be  advanced  in  the  coming  years  and  provide  increased  speed,  and  performance  for  ICs. 
However,  as  device  dimensions  continue  to  shrink,  quantum  phenomena  are  going  to  become  the 
dominant  device  factor.  There  are  currently  only  vague  ideas  on  the  use  of  quantum  wells  for  storage 
of  electrons  and  resonant  tunneling  between  wells  as  the  dominant  transport  mechanisms  between 
regions,  it  is  time  to  begin  to  seriously  investigate  such  concepts.  We  have  defined  a  research 
program  to  investigate  such  concepts  and  to  develop  analytical  techniques  to  provide  the  necessary 
characterization  and  feedback  for  understanding  structures  with  10-100  Angstrom  dimensions.  This 
program  is  outlined  in  detail  in  the  proposal  for  the  follow-on  JSEP  program. 
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In  order  to  study  quantum  well  structures,  we  must  first  understand  the  physics  and  control  of  the 
AIGaAs/GaAs  interface,  particularly,  the  asymmetries  which  are  dependent  on  the  growth  sequence. 
To  begin  developing  a  process  capability  for  lll-V  materials  and  to  investigate  the  AIGaAs/GaAs 
interface,  we  selected  a  MODFET  structure  as  a  starting  point.  This  structure  is  based  upon  the 
GaAs/AIGaAs  interface,  requires  MBE  growth  capability  and  is  very  sensitive  to  growth  and  process 
induced  changes  in  the  interface. 

A  mask  set  was  designed  with  the  aim  of  being  able  to  measure  a  whole  host  of  parameters  on  the 
same  die  with  MODFET  structures.  Chip  layout  is  such  that  the  left  side  of  the  die  can  be  split  off, 
bonded  and  cooled  to  liquid  nitrogen  or  helium  temperatures.  Of  primary  importance  on  this  side  are 
the  Hall  bars,  one  of  which  is  gated.  Also  included  are  fat  gate  and  small  geometry  MODFETs.  On  the 
other  half  of  the  die  are  more  MODFETs,  a  large  pad  for  C-V  or  DLTS  measurements,  and  an  unused 
area  for  optical  characterization,  such  as  photoluminesence.  Features  exist  on  this  side  for 
determining  various  process  parameters  such  as  implant  activation,  contact  resistance,  quality  of 
isolation  etc.  Two  thirty-nine  element  ring  oscillators,  one  for  depletion  and  one  for  enhancement 
style  MODFETs,  each  with  their  own  output  are  also  included. 

With  this  mask  set  it  should  be  possible  to  determine  the  mobility  of  the  two-dimensional  electron 
gas  at  the  heterojunction  interface  as  a  function  of  both  temperature  and  gate  bias  for  a  given  sample. 
Both  the  thickness  of  the  AIGaAs  layer  and  type  of  surface  passivation  layer,  if  any,  can  also  be  varied 
from  wafer  to  wafer. 

From  the  gated  Hall  bar  measurements  and  other  measurements  available  on  the  die,  we  hope  to 
determine  the  effect  of  filling  different  quantum  levels  at  the  interface.  This  investigation  will  include  a 
study  of  the  mobility  of  the  carriers  when  the  upper  levels  become  filled  and  the  scattering  processes 
which  are  effective  in  the  various  regimes.  A  second  effort  is  in  progress  to  predict  the  two- 
dimensional  electron  gas  mobility  using  scattering  theory  under  various  conditions  of  doping  and 
temperature. 
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2.  INFORMATION  SYSTEMS 

2.1  Real  Time  Statistical  Signal  Processing 

Principal  Investigator:  T.  Kailath 

2.1.1  Scientific  Objectives 

The  goals  of  this  research  program  were  to  study  emitter  location  problems  and  in  particular  focus 
on  methods  for  determining  the  number  of  sources  and  processing  array  data  when  the  signal  is 
completely  coherent  with  one  or  more  interferences. 

2.1 .2  Progress 

2.1 .2.1  Determining  the  Number  of  Sources 

A  new  approach  to  the  problem  of  determining  the  number  of  sources  for  a  passive  array  was 
developed  [6].  The  method  is  based  on  information  theoretic  criteria  concepts  for  model  selection 
introduced  by  Akaike  (A  Information  Criteria  -  AIC)  and  Rissanen  (Minimum  Description  Length 
•  MDL)  .The  conventional  approach  for  determining  the  number  of  sources  in  the  eigenstructure 
framework  were  based  on  hypothesis  testing  that  required  definition  of  thresholds.  The  new 
approach  does  not  require  any  subjective  threshold  settings.  The  number  of  sources  is  determined 
merely  by  minimizing  the  AIC  or  MDL  criteria. 

It  was  analytically  proven  that  the  MDL  criteria  yields  consistent  estimates  for  the  number  of 
sources  whereas  the  AIC  overestimates  their  number. 

Simulation  studies  were  carried  out  to  verify  these  methods  and  yielded  good  results  under  varying 
situations  including  low  signal  to  noise  ratios,  small  number  of  data  samples  and  closely  spaced 
signals. 

The  new  approach  is  equally  applicable  in  other  areas  where  eigenstructure  methods  have  been 
used.  Some  examples  are  determining  the  number  sinusoids  in  time  series  data,  the  number  of  poles 
from  the  natural  response  of  a  linear  system  and  the  number  of  overlapped  echoes  in  backscatter 
data  typical  in  radar,  etc. 


JSEP  Annual  Report 


Page  73 


June  1984 


Stanford  University 


Stanford  Electronics  Laboratories 


2. 1.2. 2  Coherent  Signals  and  Interference 

The  key  concept  here  has  been  the  use  of  spatial  smoothing  to  decorrelate  coherent  arrivals  at  the 
sensor  array  thus  effectively  solving  the  coherent  signal  and  interference  problem.  It  is  of  interest  to 
note  that  the  problem  of  dealing  with  fully  coherent  signals  and  interference  has  been  outstanding  for 
nearly  two  decades  after  it  was  found  that  optimum  beamforming  techniques  do  not  work  in  these 
conditions. 

In  [5]  we  show  how  spatial  smoothing  can  be  used  to  restore  the  rank  of  the  covariance  matrix  of 
the  signals  and  coherent  interference  impinging  on  an  array.  This  enables  the  use  of  eigenstructure 
methods  to  determine  the  directions  of  arrival  of  sources  while  still  yielding  estimates  which  are 
unbiased  and  asymptotically  exact. 

In  [3]  and  [4]  we  applied  these  ideas  to  beamforming  and  showed  how  to  carry  out  optimum 
beamforming  in  a  coherent  signal  environment  without  suffering  signal  cancellation.  Simulation 
studies  were  carried  out  for  various  scenarios  including  narrowband  and  broadband  sources, 
different  signal  to  noise  ratios  and  several  source  and  sensor  directions  to  assess  the  performance  of 
the  new  beamformer.  The  results  have  borne  out  theoretical  predictions.  Signal  cancellation  has  been 
eliminated  and  effective  nulling  of  coherent  jammers  obtained. 

Some  attractive  features  of  spatial  smoothing  concept  include  simplicity,  low  computational 
overhead  and  robustness  to  actual  degree  of  coherence. 

2.1 .2.3  Other  Applications  of  Eigenstructure  Methods 

The  eigenstructure  methods  have  been  applied  to  new  problems  in  statistical  signal  processing.  In 
[1]  we  discussed  new  methods  to  resolve  highly  overlapping  echoes.  The  signal  model  is  typical  in 
radar  and  geophysics  where  repetitive  observations  are  available  of  overlapped  signals  in  additive 
noise.  The  performance  of  this  new  method  has  been  most  promising.  Applications  of  this  high 
temporal  resolution  technique  to  radar/sonar  target  imaging  (for  identification)  and  medical  imaging 
are  under  investigation. 

Another  study  [7]  addresses  source  localization  by  widely  spaced  subarrays.  Processing 
techniques  which  operate  only  on  subarray  covariances  have  been  developed  which  outperform 
conventional  triangularization  methods.  Moreover,  consistent  estimates  for  the  number  of  sources  are 
obtained. 

In  [2],  some  new  concepts  were  reported  for  optimal  beamforming  in  situations  where  one  type  of 
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interference  (say  directional  jammers)  are  considered  more  undesired  than  say  sensor  to  sensor 
uncorrelated  noise.  Eigenstructure  methods  have  been  used  to  build  a  modified  beamformer  with 
these  desired  properties. 


2.1.3  Summary  of  Results 

1.  Derivation  of  a  consistent  information  theoretic  estimator  for  a  number  of  sources. 

2.  A  new  method  to  decorrelate  coherent  wavefronts  which  can  be  applied  to  a  wide  range 
array  processing  methods  to  effectively  solve  this  long  outstanding  problem. 

3.  New  applications  of  eigenstructure  concepts  to  resolving  overlapping  signals,  source 
localization  using  subarray  covariances  and  optimum  processing  when  some 
interferences  are  more  undesired  than  others. 
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2.2  Signal  Processing  and  Compression 

Principal  Investigators:  R.M.  Gray,  M.  E.  Heilman,  T.M.  Cover,  J.  Gill 
2.2.1  Scientific  Objectives 

The  objective  of  this  project  is  to  study  the  basic  limitations  on  the  performance  and  efficiency  of 
communication  networks,  including  both  multiuser  and  point-to-point  systems.  A  parallel  objective  is 
to  develop  improved  methods  of  signal  processing  for  these  systems  from  the  structures  used  to 
verify  the  performance  limitations.  Attention  will  be  directed  toward  the  trade-offs  between 
performance  and  complexity  in  signal- processing  algorithms.  Signal  processing  for  data 
compression,  improved  signal-to-noise  ratios,  and  classification  will  be  investigated. 


2.2.2  Progress 

2.2.2. 1  Vector  Quantization 
investigators:  (R.M.  Gray,  P.C.  Chang,  T.  Flynn) 


During  the  past  year  the  project  has  focused  on  three  principal  problems:  The  design  of  feedback 
vector  quantizers  for  data  compression,  the  design  of  adaptive  vector  quantizers  for  speech 
compression,  and  the  extension  of  vector  quantizer  design  algorithms  to  the  design  of  combined 
quantization/classification  distributed  sensor  networks. 

•  Feedback  Vector  Quantizers: 

A  stochastic  gradient  algorithm  for  designing  predictive  vector  quantizers  (PVQ)  was 
developed  in  order  to  get  better  codebooks  and  faster  design  speed.  To  further  improve 
the  performance,  adaptive  VQ's  which  combine  the  waveform  coding  techniques  and 
linear  predictive  coding  (LPC)  techniques  were  studied. 

Predictive  Vector  Quantizers  are  vector  generalizations  of  scalar  predictive  quantization 
or  DPCM:  a  vector  quantizer  in  a  feedback  loop  quantizes  the  error  between  the  new 
input  vector  and  a  reproduction  vector  predicted  from  past  quantizer  outputs.  Both  the 
predictor  and  quantizer  operate  on  vectors.  Because  of  the  feedback  loop,  it  is  more 
complicated  to  design  a  locally  optimum  codebook  for  a  PVQ  than  for  a  memoryless  VQ. 
Cuperman  and  Gersho  [3, 2,  5]  designed  a  PVQ  with  two  main  steps:  first  a  set  of  linear 
vector  predictive  coefficients  was  designed  by  generalized  LPC  techniques;  then  a 
codebook  for  these  predictor  coefficients  was  designed  by  the  generalized  Lloyd 
algorithm.  The  codebook  obtained  by  this  method  may  be  locally  optimum  under  the 
assumption  of  very  small  quantizer  error.  Unfortunately,  this  assumption  is  not  valid  for 
low  rate  coding.  Furthermore,  no  procedures  were  given  to  improve  the  predictor  for  the 
quantizer,  therefore,  once  the  predictor  coefficients  are  not  optimum,  the  final  codebook 
is  nonoptimum. 
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Codebook  and  predictor  design  based  on  a  stochastic  gradient  algorithm  (or, 
equivalently,  a  generalized  steepest  descent  algorithm)  avoids  these  disadvantages:  It 
simultaneously  improves  the  linear  predictive  coefficients  and  the  vector  quantizer  in  the 
feedback  loop  iteratively  for  a  given  training  sequence.  This  provides  better  performance 
and  faster  convergence.  Preliminary  results  may  be  found  in  Gray  [6]. 

•  Adaptive  Vector  Quantization  for  Speech: 

An  adaptive  vector  quantizer  is  a  system  using  a  model  classifier  to  adapt  a  waveform  VQ 
so  that  the  waveform  codebook  can  be  better  matched  to  local  behavior  of  speech, 
thereby  achieving  better  reproduction  quality.  Both  the  adaptation  information  and 
waveform  encoding  result  are  sent  to  the  receiver.  The  model  VQ  typically  operates  on  a 
much  larger  vector  of  samples,  e.g.,  128  samples,  and  the  codebook  size  of  the  model  VQ 
is  small,  hence  the  number  of  bits  spent  on  specifying  the  model  through  the  side  channel 
are  typically  much  less  than  those  devoted  to  the  waveform  coder. 

We  developed  a  simulated  code  design  algorithms  for  several  adaptive  VQ's  using  a  voice 
coding  LPC  VQ  with  an  Itakura-Saito  distortion  and  ordinary  memoryless  VQ’s  with 
squared  error  distortion  for  the  data  classified  by  the  LPC  VQ.  Preliminary  results  and 
comparisons  with  other  compression  systems  may  be  found  in  Gray  [6]. 

•  Quantization  and  Classification  in  Distributed  Sensor  Networks: 

This  work  treats  the  problem  of  data  quantization  in  distributed  sensor  (e.g.,  radar  or 
sonar)  systems.  The  sensors,  at  several  widely  separated  locations,  gather  data  which 
must  be  converted  to  digital  form  for  communication  to  processing  locations  at  which  the 
digitized  data  are  combined  and  the  inference  (e.g.,  detection,  location,  or  classification) 
is  done.  The  digitized  data  must  be  as  informative  as  possible,  subject  to  the  limitations 
on  data  rate.  The  appropriate  strategy  is  not  to  reproduce  the  original  data  waveform,  but 
to  preserve  its  sensitivity  to  the  targets  under  observation. 

The  problems  of  detection  and  classification  (in  which  one  chooses  among  a  finite 
number  of  alternatives)  were  formulated  in  a  Bayesian  fashion.  Two  methods  of  quantizer 
design  were  considered:  one  which  makes  small  the  exact  expected  cost,  the  other 
which  makes  small  an  upper  bound  expressed  as  a  weighted  sum  of  distributional  overlap 
measures.  The  second  method  was  motivated  by  the  guidelines  it  provides  for 
cooperation  among  sensor  locations.  In  both  methods,  a  quantizer  is  specified  by  a 
"codebook"  of  representative  points  and  a  distortion  measure  which  dictates  how  data 
points  are  mapped  into  codewords;  this  form  allows  a  prototype  quantizer  to  be  improved 
using  an  iterative  algorithm  similar  to  the  Lloyd  algorithm. 

In  problems  of  estimation  (of  a  real  vector  parameter),  attention  was  restricted  to  cases 
where  the  log-likelihood  function  induced  by  an  observation  can  be  treated  as  quadratic. 
(This  is  true,  for  example,  for  a  strong  known  signal  contaminated  with  Gaussian  noise.) 
described  in  terms  of  them.  The  increase  in  mean-square  estimation  error  due  to 
quantization  was  found  to  be  expressible  (approximately)  as  a  weighted  mean-square 
error  in  the  quadratic  parameters.  Using  this,  the  established  methods  of  quantization  for 
minimum  mean-square  error  can  be  applied  to  this  problem. 
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2.2.2.2  Advances  in  Cryptography 

Investigator:  M.E.  Heilman 

In  [4]  we  presented  a  public  key  distribution  system  which  depended  on  the  difficulty  of  computing 
discrete  logarithms.  Merkle  [7]  and  [1]  have  found  a  faster,  subexponential  method  for  computing 
this  function.  In  earlier  work,  we  extended  this  algorithm  from  GF(q)  to  the  more  general  case, 
GF(q**m)  with  q  fixed  and  m  tending  to  infinity.  In  doing  this,  we  also  found  a  way  to  rigorize  the  use 
of  a  precomputation,  by  defining  a  "virtual  spanning  set". 

In  recent  work,  we  have  extended  the  subexponential  behavior  to  GF(q*  *2).  This  was  an  important 
step  because  it  indicates  that  the  discrete  logarithm  function  is  probably  computable  in 
subexponential  time  in  general.  At  first  we  were  only  able  to  show  the  subexponential  behavior  for 
99.8%  of  the  primes  q  because  not  ail  quadratic  fields  possess  the  unique  factorization  property.  But, 
by  defining  "super  primes"  to  be  quadratic  primes  which  also  have  prime  norm,  we  were  able  to 
circumvent  this  problem  and  extend  the  result  to  all  primes  q. 

We  also  have  developed  a  new  digital  signature  scheme  based  on  the  discrete  logarithm  problem. 
The  signature  scheme  differs  from  the  RSA  system  in  that  it  is  not  based  on  a  public  key 
cryptosystem. 

2. 2. 2. 3  Advances  in  Data  Compression  and  Information  Inequalities 

Investigator:  Thomas  M.  Cover 


Our  work  covers  several  areas. 

•  Multiple  Descriptions: 

In  Cover  and  El  Gamal,  we  present  a  rate  region  for  the  multiple  descriptions  problem. 
Here  we  wish  to  give  two  separate  descriptions  of  a  random  variable  such  that  each  is 
good  but  together,  they  are  very  good.  Recently  we  have  obtained  partial  results  on  the 
converse  for  such  a  rate  region  with  Shen  Shi  Yi.  We  await  further  developments.  Berger 
and  Zhan  Zheng  have  informally  announced  some  results  showing  that  the  Cover/El 
Gamal  region  is  only  an  inner  bound  rather  than  the  true  rate  distortion  rate  region  for 
this  problem. 

•  Hadamard: 

We  have  found  an  unsuspected  new  proof  of  the  classical  Hadamard  inequality  that  says 
that  the  determinant  of  a  matrix  is  less  than  the  products  of  the  lengths  of  its  rows.  The 
three  line  proof  follows  from  plugging  the  multivariate  normal  into  the  entropy  inequality 
and  represents  the  first  new  proof  of  this  theorem  in  30  years. 
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•  AEP 

Our  work  on  showing  that  there  is  an  asymptotic  equipartition  property  for  multiplicative 
processes  (like  the  stock  market)  is  still  in  progress.  The  AEP  first  became  known  in 
statistical  mechanics  and  was  later  used  to  good  effect  by  Shannon  in  information  theory. 
It  states  that  with  high  probability,  an  ergodic  source  really  consists  of  an  almost  uniform 
distribution  over  a  small  number  of  possible  outcomes,  hence  the  term  "equipartition." 
Nicely  enough,  such  an  equipartition  characterization  is  also  true  of  an  ergodic  stock 
market.  It  is  shown  that  the  capital  growth  rate  (rather  than  the  entropy  rate)  tends  to  a 
limit  characteristic  of  the  market.  This  allows  comparing  compounding  random 
processes  by  their  growth  rates  rather  than  by  looking  deeply  into  the  guts  of  the 
statistical  time  dependence. 

•  Brunn  Minkowski: 

In  the  paper  Costa/Cover,  we  show  a  direct  analogy  between  the  Brunn-Minkowski 
inequality  for  the  volumes  of  set  sums  under  addition  and  the  deep  entropy  power 
inequality.  The  similarity  is  so  striking  that  there  must  exist  a  common  proof.  Although 
we  do  not  find  such  a  proof,  we  exhibit  new  theoretical  consequences  of  both 
inequalitities. 

•  Burg's  Theorem: 

In  Choi/Cover,  currently  scheduled  to  appear  in  IEEE  Proceedings,  we  find  a  new  four 
line  proof  that  Burg's  maximum  entropy  spectral  density  estimate  is  achieved  by  the 
appropriate  Gauss- Markov  process.  This  simplifies  previous  proofs  which  used  irrelevant 
mathematics  to  establish  the  result. .  The  proof  is  primarily  information-theoretic  and 
follows  from  the  chain  rule  for  entropies. 


2. 2. 2.4  Data  Compression  for  Computer  Data  Structures 
investigators:  J.  Gill,  K.  Chu 


In  the  past  year,  this  project  has  evolved  from  primarily  theoretical  and  software  efforts  concerned 
with  compacting  tree-like  data  structures.  Our  research  has  focused  on  the  development  of  a 
hardware/software  subsystem  for  compression  of  computer  file  systems.  The  goal  of  the  subsystem, 
whose  preliminary  design  was  completed  in  February  1984,  is  to  achieve  a  compression  ratio  of  50% 
at  a  data  rate  of  1  megabit/second.  Although  the  specified  data  rate  is  an  order  of  magnitude  slower 
than  that  of  typical  computer  disk  subsystems,  it  is  two  orders  of  magnitude  faster  than  that 
achievable  by  software-only  approaches. 


The  data  compression  technique  we  have  selected  is  a  modification  of  the  method  of  Ziv  and 
Lempef  [8, 9J.  The  implementation  is  an  improvement  of  the  suffix-tree  approach  of  Rodeh,  Even,  and 
Pratt.  The  improved  algorithm  yields  faster  execution  for  the  micro-coded  compression  engine  that 
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we  have  designed,  although  the  resulting  compression  ratio  is  somewhat  poorer.  The  compression 
and  expansion  algorithms  were  first  implemented  in  software  using  C  programs  running  on  the  Vax 
11/780  in  the  Information  Systems  Laboratory.  In  tests  on  more  than  5  megabytes  of  text  and  binary 
data,  we  confirmed  that  a  compression  ratio  of  slightly  more  than  50%  can  be  achieved. 

An  evaluation  of  the  software  versions  of  these  compression  algorithms  revealed  that  the  great 
majority  of  the  execution  time  was  spent  in  die  relatively  simple  operation  of  following  links  in  a  suffix 
tree.  Since  these  operations  can  be  performed  very  rapidly  by  low  cost  special  purpose  hardware,  we 
began  a  study  of  the  feasibility  of  a  hardware  implementation  of  our  improved  suffix-tree  construction. 
The  basic  hardware  design,  based  on  the  AMO  29116  16-bit  bipolar  microcontroller  and  the  AMD 
2910  microprogram  sequencer,  was  completed  in  February,  as  was  the  micro-instruction  format  for 
the  compression  engine.  The  prototype  design  comprises  approximately  70  integrated  circuits;  most 
of  the  ICs  are  memory  and  latches,  which  allow  the  microprogram  to  be  stored  in  RAM  during 
subsystem  development.  The  prototype  compression  unit  is  standalone;  a  revised  version  will  be 
attached  to  the  data  bus  of  a  small  computer  (IBM  PC  compatible),  so  that  data  transfers  to/from  the 
disk  drive  can  be  routed  through  the  compression  unit. 

In  summary,  algorithms  for  compression  of  data  in  computer  file  systems  have  been  developed  and 
implemented  in  software.  The  preliminary  design  of  a  low-cost  hardware  implementation  has  been 
completed. 

3.  Summary  of  Results 

Two  new  code  structures  and  design  techniques  for  these  strucures  were  developed  for  data 
compression  systems:  predictive  vector  quantizers  designed  by  stochastic  gradient  algorithms  and 
adaptive  vector  quantizers  using  memoryless  vector  quantizers  for  both  the  classifier  used  to  adapt 
and  for  the  waveform  coder.  A  new  approach  was  developed  to  quantizing  data  in  distributed  sensor 
networks  so  as  to  best  communicate  the  information  relevant  to  other  nodes'  required  decisions. 
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